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ABSTRACT 

The major uncertainties involved in the Chandrasekhar mass models for Type la supernovae (SNe 
la) are related to the companion star of their accreting white dwarf progenitor (which determines the 
accretion rate and consequently the carbon ignition density) and the flame speed after the carbon ignition. 
We calculate explosive nucleosynthesis in relatively slow deflagrations with a variety of deflagration speeds 
and ignition densities to put new constraints on the above key quantities. The abundance of the Fe- 
group, in particular of neutron-rich species like ''^Ca, ^°Ti, ^^Cr, ^"''^^Fe, and ^^Ni, is highly sensitive to 
the electron captures taking place in the central layers. The yields obtained from such a slow central 
deflagration, and from a fast deflagration or delayed detonation in the outer layers, are combined and 
put to comparison with solar isotopic abundances. To avoid excessively large ratios of ^"^Cr/^^Fe and 
^°Ti/^^Fe, the central density of the "average" white dwarf progenitor at ignition should be as low as 2 
X 10^ g cm"'^. To avoid the overproduction of ^*Ni and ^"^Fe, either the flame speed should not exceed 
a few % of the sound speed in the central low Yg layers, or the metallicity of the average progenitors 
has to be lower than solar. Such low central densities can be realized by a rapid accretion as fast as 
M ^1 y. 10~''Mq yr~^. In order to reproduce the solar abundance of ''^Ca, one also needs progenitor 
systems that undergo ignition at higher densities. Even the smallest laminar flame speeds after the 
low-density ignitions would not produce sufficient amount of this isotope. We also found that the total 
amount of ^^Ni, the Si-Ca/Fe ratio, and the abundance of some elements like Mn and Cr (originating 
from incomplete Si-burning) , depend on the density of the deflagration-detonation transition in delayed 
detonations. Our nucleosynthesis results favor transition densities slightly below 2.2x10^ g cm^'^. 



1. INTRODUCTION 

There are strong observational and theoretical indica- 
tions that Type la supernovae (SNe la) are thermonuclear 
explosions of accreting white dwarfs (e.g., Wheeler et al. 
1995; Nomoto, Iwamoto & Kishimoto 1997; Branch 1998). 
Theoretically, both the Chandrasekhar mass white dwarf 
models and sub-Chandrasekhar mass models have been 
considered (see, e.g., Arnett 1996; Nomoto et al. 1994, 
1996a, 1997b, 1997c; Canal, Ruiz-Lapuente, & Isern 1997 
for reviews of recent progress) . Though these white dwarf 
models can account for the basic observational features 
of SNe la, the exact binary evolution that leads to SNe la 
has not been identified yet. Various evolutionary scenarios 
have been proposed, which include (1) a double degenerate 
scenario, i.e., the merging of two C-f-0 white dwarfs in a 
binary system with a combined mass exceeding the Chan- 
drasekhar mass limit (e.g., Iben & Tutukov 1984; Webbink 
1984) and (2) a single degenerate scenario, i.e., accretion of 
hydrogen or helium via mass transfer from a binary com- 
panion at a relatively high rate (e.g., Nomoto 1982a). In 
the case of helium accretion at low rates, He detonates at 
the base of the accreted layer before the system reaches the 
Chandrasekhar mass (Nomoto 1982b; Woosley & Weaver 
1986, 1994a; Livne & Arnett 1995). Currently, the is- 
sues of the Chandrasekhar mass versus sub-Chandrasekhar 



mass models and the double degenerate versus single de- 
generate scenarios are still debated (see, e.g., Renzini 1996 
and Branch et al. 1995 for recent reviews), but they are 
being confronted with an increasing number of observa- 
tional constraints. 

The observational search for the double degenerate sce- 
nario led to the discovery of a few binary white dwarfs 
systems, but with combined mass being smaller than the 
Chandrasekhar mass (Renzini 1996). Theoretically, it has 
been suggested that the merging of double white dwarf 
systems leads to accretion-induced collapse rather than to 
SNe la (Nomoto & Iben 1985; Saio & Nomoto 1985, 1998). 
The Chandrasekhar versus sub-Chandrasekhar mass issue 
has recently experienced some progress. Photometric and 
spectroscopic features of SNe la in early phases clearly 
indicate that Chandrasekhar mass models give a much 
more consistent picture than the sub-Chandrasekhar mass 
models of helium detonations (e.g., Hoflich & Khokhlov 
1996; Nugent et al. 1997). This leaves us with the most 
likely progenitor system, a single degenerate system with 
hydrogen accretion from the companion star, leading to 
a Chandrasekhar-mass white dwarf. However, the Chan- 
drasekhar mass model W7 (Nomoto, Thielemann, & Yokoi 
1984; Thielemann, Nomoto, & Yokoi 1986), widely used in 
galactic chemical evolution calculations, may require im- 
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provements in terms of the Fc-group composition because 
it predicts significantly higher ''''^Ni/'''^'Fe ratios tfian solar. 
The direct determination of Ni abundances in late time 
SN la spectra is therefore important (Ruiz-Lapuente 1997; 
Liu, Jcffcry, & Schultz 1997; Mazzah et al. 1998). The 
recent findings of supersoft X-ray sources, being potential 
progenitors of SN la events with high accrc;tion rates, caus- 
ing ignition at low densities (van den Heuvel et al. 1992; 
Rappaport, Di Stefano, fc Smith 1994; Di Stefano et al. 
1997), leave hope for Chandrasckhar-mass models, which 
meet all these requirements, also for late time spectra. 

The presupernova evolution of an accreting white dwarf 
depends on the accretion rate M, the composition of the 
material transferred from the companion star, and the ini- 
tial mass of the white dwarf (e.g., Nomoto 1982a; Nomoto 
& Kondo 1991). Chandrasekhar mass white dwarfs can 
be obtained with a relatively high mass transfer rate of 
hydrogen of the order M Ki A x 10~^ - 10~^ M© yr~^. 
At M > 4 X 10~^ Mq yr~^, the accreting white dwarf 
blows off a strong wind, which reduces M to an effective 
accretion rate below lO^*' Mq yr^^ (Hachisu, Kato, & 
Nomoto 1996, 1999a). This avoids the formation of an 
extended envelope in the accreting white dwarf. (Nomoto 
, Nariai, & Sugimoto 1979). At such rates hydrogen and 
helium burn steadily or with weak flashes, leading to a 
white dwarf with a growing C+0 mass. 

For the Chandrasekhar mass white dwarf model, carbon 
ignition in the central region leads to a thermonuclear run- 
way. The ignition density depends on the stellar structure 
as a function of previous accretion history. High accretion 
rates lead to higher central temperatures, i.e. favoring 
lower ignition densities. A flame front then propagates 
at a subsonic speed v^ef as a deflagration wave owing to 
heat transport across the front (Nomoto et al. 1984). The 
major and yet not fully solved questions are related to 
the propagation of the burning front. Timmes & Woosley 
(1992) have analyzed the propagation speed of laminar 
burning fronts in one dimension as a function of density 
and fuel composition. However, the propagation in three 
dimensions is influenced by instabilities that can enhance 
the effective radial flame speed beyond its laminar value. 

The flame front is subject to various types of insta- 
bilities, namely, thermal instabilities (Bychkov &: Liber- 
man 1995a), the Landau-Darrius (L-D) instability (Lan- 
dau & Lifshitz 1987), the Rayleigh-Taylor (R-T) instabil- 
ity, and the Kevin- Helmholtz (K-H) instability (Niemeyer, 
Woosley, & Hillebrandt 1996). The turbulent burning 
regime associated with the R-T bubbles on global scales 
has been studied (Livne 1993; Arnett & Livne 1994a; 
Khokhlov 1995; Niemeyer & Hillebrandt 1995a; Niemeyer 
et al. 1996), but there remain many uncertainties, re- 
lated partially to numerical resolution but also to the role 
and spectrum of turbulent length scales (Hillebrandt & 
Niemeyer 1997). With the present uncertainties, it is es- 
sential to perform parameterized sets of calculations that 
explore the possible range of effective radial flame speeds. 

In the deflagration wave, electron captures enhance the 
neutron excess. The amount of electron capture depends 
on both v^ei (influencing the time duration of matter at 
high temperatures, and with it the availability of free pro- 
tons for electron capture and the high-energy tail of the 
electron energy distribution) and the central density of 



the white dwarf = pc/lO^ g cm^^ (increasing the elec- 
tron c;hcmical potential). The resultant nucleosynthesis in 
slow deflagrations (see, e.g., Khokhlov 1991b) has some 
distinct features compared with faster deflagrations like 
W7 (Nomoto et al. 1984; Thielemann et al. 1986), thus 
providing important constraints on these two parameters. 
The constraint on the central density is equivalent to a 
constraint on the accretion rate, as discussed above. After 
an initial deflagration in the central layers, the deflagra- 
tion is accelerated as in W7, or assumed to turn into a 
detonation at lower densities, as in the delayed detona- 
tion models (Khokhlov 1991a; Woosley & Weaver 1994a). 
For the latter, the uncertain transition density ptr would 
result in a variety of total masses of ^^Ni and expansion 
velocities of the outer layers. 

To obtain constraints on the three parameters (pigrn 
Wdef) and ptr), we performed explosive nucleosynthesis cal- 
culations for slow deflagrations followed by a delayed det- 
onation or a fast deflagration. These calculations assumed 
spherical symmetry and therefore might not be fully ade- 
quate for a realistic and consistent approach, but we ex- 
pect that they lead at least to some clues how abundance 
features relate to the spherical average of these quantities 
in realistic models. Initially slower deflagrations cause an 
earlier expansion of the outer layers with respect to the 
arrival of the burning front (as information of the cen- 
tral ignition propagates with sound speed; Nomoto, Sugi- 
moto, & Neo 1976) and lead to low densities for the outer 
deflagration and detonation layers. In such a case, even 
a detonation does not lead to a pure Fe-group composi- 
tion, as expected in central detonations, and intermediate 
mass elements from Si to Ca are produced at a wide range 
of expansion velocities. Therefore, if the deflagration- 
detonation transition (DDT) density is well tuned, delayed 
detonations can meet these observational requirements as 
well as fast deflagration (Kirshner et al. 1993). Compared 
with the earlier delayed detonation models by Khokhlov 
(1991b), we adopt a larger and more detailed nuclear re- 
action network that alos includes electron screening. In 
comparison to the fast deflagration models by Woosley 
(1997b), we use initial models with lower central densities 
and smaller flame speeds, to concentrate on our main aim, 
which is to find the " average" SN la conditions responsible 
for their nucleosynthesis contribution to galactic evolution, 
i.e., especially the Fe-group composition. 

From the very early days of explosive nucleosynthesis 
calculations, when no direct connection to astrophysical 
sites was possible yet, it was noticed (Trimble 1975) that 
the solar Fe-group composition could be reproduced with 
a superposition of matter from explosive Si burning with 
about 90% originating from a Ye = 0.499 source and 10% 
from a Ye = 0.46 source. As it has been shown that SN 
II ejecta with Yg < 0.498 could cause serious problems 
in comparison with observations (Thielemann, Nomoto, & 
Hashimoto 1996), SNe la have to be identified with this 
second source and the appropriate conditions which lead 
to the best agreement with solar abundances. For this rea- 
son we present detailed yields of delayed detonation mod- 
els as well as fast deflagration models and compare them 
with solar abundances for a number of "training sets" of 
ignition densities, flame speeds, and DDT densities. 

These nucleosynthesis constraints can provide clues to 
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the explosion mechanism (i.e., the speed of the burning 
front) and the ignition density (i.e., the accretion rate fr'om 
the binary companion) for the " average" or dominating SN 
la contributions during galactic chemical evolution. This 
is the purpose of the present paper. We have to be aware 
that there are some systematic variations that manifest 
themselves in light curves (i.e., the brighter one is slower; 
Phillips et al. 1990; Hamuy et al. 1995) and might lead to 
a variation in nucleosynthesis as well (Hoflich & Khokhlov 
1996, Hoflich, Wheeler & Thielemann 1998). It is further 
of importance to explore metallicity effects, which might 
have an influence on the evolution of the progenitor sys- 
tems, with respect to the initial mass function (IMF) and 
composition of white dwarfs as well as the binary accretion 
history (e.g., Yoshii, Tsujimoto, & Nomoto 1996; Hoflich 
et al. 1998; Umeda et al. 1999; Hoflich et al. 1999). 
Only the latter will clarify whether the nature of SNe la 
at high redshifts is the same as for nearby SNe la, which 
enters the determination of the cosmological parameters 
Ho and qo (e.g.. Branch & Tammann 1992; Riess, Press, 
& Kirshner 1995; Riess et al. 1999; Perlmutter et al. 1997, 
1999). Observable spectral features could possibly help to 
identify the metallicity due to slight changes in nucleosyn- 
thesis (Hoflich et al. 1998; Hatano et al. 1999; Lentz et 
al. 1999). 

After a description of our model calculations in §2, in- 
cluding initial models, the hydrodynamic treatment, and 
a discussion of the input physics, we present in §3 detailed 
nucleosynthesis results from slow deflagrations and delayed 
detonations in comparison with the carbon deflagration 
model W7. In §4, the integrated abundances of SNe la 
models are combined with those of SNe II to compare with 
solar abundances. Finally, we discuss constraints on possi- 
ble evolutionary scenarios and give conclusions in §5. Very 
preliminary accounts of the present investigations on nu- 
cleosynthesis in slow deflagrations and delayed detonations 
have been given in Nomoto et al. (1997c) and Thielemann 
et al. (1997). 

2. INITIAL MODELS, EXPLOSION HYDRODYNAMICS, AND 
INPUT PHYSICS 

2.1. Initial Models 

We adopt two models with central densities of pg = 1.37 
(C) and 2.12 (W) at the onset of thermonuclear runaway, 
i.e., at the stage when the timescale of the temperature 
rise in the center becomes shorter than the dynamical 
timescale. Here C and W imply that these are the same 
models as calculated for C6 and W7, respectively (Nomoto 
et al. 1984). The initial white dwarfs of these models, be- 
fore the onset of H-accretion, have a mass of M = 1.0 Mq, 
a central temperature of Tc = 1.0 x lO'^ K, and com- 
positions of XC^G) = 0.475, X(i60) = 0.50, X{^^Ne) = 
0.025. 

The outer layers of the mass grid extend to the steady 
hydrogen burning shell as an outer boundary. The temper- 
ature and density at the burning shell are determined from 
the boundary condition that the accreted matter is pro- 
cessed into helium with the mass accretion rate Al. These 
values increase from initally 8 x 10^ K and 1 x 10'' g 
cm~^ to 1 X 10^ K and 1 x 10® g cm^'^ at the point 
of central carbon ignition. The accretion rate for case C 
is due to steady and stable hydrogen burning correspond- 



ing to the C-f-O core increase during an asymptotic giant 
branch evolution (see Nomoto 1982a). The rate for case W 
is kept constant up to the point of carbon ignition at the 
center. The exact values are given by equations (1) and 
(2), where M indicates the mass of the accreting white 
dwarf. 

M(C) = 8.5 X lO"^(Af/M0 - 0.52)MQy"i (1) 
M(W) = 4 x lO^^^Mgy-^ (2) 

During the accretion phase the white dwarf mass M in- 
creases with time and the central temperature increases as 
a result of heat inflow from the H-burning shell as well as 
compressional heating. Cooling is due mostly to plasmon 
neutrino losses and neutrino bremsstrahlung. Cooling due 
to Urea shells and the convective Urea process is not taken 
into account. This has no effect on case C, but could delay 
the ignition in case W to higher densities (e.g., Paczynski 
1973; Iben 1982; Nomoto & Iben 1985; Barkat & Wheeler 
1990). When the central density pc reaches 1.5 x 10^ g 
cm"-^ (C) or 2.5 x 10^ g cm~'^ (W), carbon is ignited in 
the center, where the nuclear energy generation rate ex- 
ceeds the neutrino losses. When the central temperature 
increases owing to carbon ignition, a convective core de- 
velops. The convective energy transport is calculated in 
the framework of the time-dependent mixing length the- 
ory (Unno 1968). At Tc 8 x 10^ K convection can no 
longer transport energy in our model and the central re- 
gion undergoes a thermonuclear runaway with pc = 1.37 
X 10^ gcm-3 (C) or 2.12 x 10^ g cm'^ (W). 

2.2. Slow Deflagrations 

We know the absolute lower limit for the deflagration 
speed after central ignition from the one-dimensional anal- 
ysis of laminar flame fronts (Timmes & Woosley 1992), 
being close to 1% of the local sound speed Vg. On the 
other hand, any hydrodynamic instability can enhance 
the speed. Our parametrized "fast" deflagration studies, 
which reached 10%-30% of the sound speed and produced 
the model W7 (Nomoto et al. 1984; Thielemann et al. 
1986), however, resulted in problematic Fe-group nucle- 
osynthesis (see also discussion below). The flame speed 
found in multidimensional hydrodynamic simulations is 
still subject to large uncertainties, ranging from fdcf/^s 
0.015 (Niemeyer & Hfllebrandt 1995a) to Vdd/vs ^ O.l 
(Niemeyer et al. 1996). Therefore, it is important to in- 
vestigate how the nucleosynthesis outcome depends on the 
flame speed. In order to contrast our "fast" deflagration 
model W7, we study "slow" deflagrations here and choose 
the following parameter ranges: After the central ther- 
monuclear runaway, we assume that a slow (S) deflagration 
propagates with speeds Vdef/vs — 0.015 (WS15, CS15), 
0.03 (WS30, CS30), and 0.05 (CS50) and consider also 
extreme cases of fully and initially laminar flame fronts 
(WLAM, WSL). The location of the deflagration wave in 
radial mass coordinate and the changes in tempera- 
ture and density are shown in Figure |l| as a function of 
time. Behind the deflagration wave, the temperature rises 
quickly to values as high as T = 9 x 10^ K and the mate- 
rial experiences nuclear statistical equilibrium (NSE). As 
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the flame front propagates outward, the white dwarf ex- 
pands slowly which reduces the central density p^. The 
decrease in pc for these slow deflagrations is significantly 
slower than in W7. 

2.3. Transition from Deflagration to Detonation 

If the deflagration speed continues to be much slower 
than in W7, the white dwarf undergoes a large amplitude 
pulsation, as first found by Nomoto et al. (1976). In this 
pulsating deflagration model, the white dwarf expands and 
nuclear burning is quenched when the total energy of the 
star is still negative. In the following contraction more ma- 
terial burns, resulting in a positive total energy E. Even- 
tually the white dwarf is completely disrupted. The model 
by Nomoto et al. (1976) resulted in i? = 5 x 10''^ ergs 
and a ^^Ni mass of Mni ~ 0.15 Mq. Such a pulsating 
deflagration produces explosion energies too small to ac- 
count for typical SNe la but might be responsible for rare 
events such as SN 1991bg. 

In order to produce sufficient amounts of radioactive 
^^Ni 0.6 Mq) to power SNe la light curves by a defiagra- 
tion wave, the flame speed must be accelerated. The de- 
gree to which the flame speed is increased depends on the 
effect of R-T instabilities during the pulsation (Woosley 
1997a). The deflagration might induce a detonation when 
reaching the low-density layers. In the delayed detona- 
tion model (Khokhlov 1991a; Woosley & Weaver 1994b), 
the deflagration wave is assumed to be transformed into 
a detonation at a specific density during the first expan- 
sion phase. In the pulsating delayed detonation model 
(Khokhlov 1991b), the transition into a detonation is as- 
sumed to occur close to the maximum compression after 
recontraction, as a result of mixing. 

Physical mechanisms by which such deflagration-to- 
detonation transitions (DDTs) occur have been studied 
by Arnett & Livne (1994b), Niemeyer & Woosley (1997), 
Khokhlov, Oran, & Wheeler (1997), and Niemeyer & Ker- 
stein (1997): (1) When a sufficiently shallow temperature 
gradient is formed in the fuel, a deflagration propagates 
as a result of successive spontaneous ignitions. Such an 
over-driven deflagration propagates supersonically. (2) If 
a sufficiently large amount of fuel has such a shallow tem- 
perature gradient, the deflagration may induce a deto- 
nation wave. The critical masses for the formation of a 
detonation are quite sensitive to the carbon mass frac- 
tion X{C), e.g., ~ 10-19 Mq and - lO^i* Mq at 
p ^ 3x 10^ g cm~'^ for X{C)=1.0 or 0.5, respectively. 
(Niemeyer & Woosley 1997). (3) Such a shallow tempera- 
ture gradient region in the fuel may be formed if the fuel 
is efficiently heated by turbulent mixing with ashes. Such 
a mixing and heat exchange may occur when the turbu- 
lent velocity associated with the flame destroys the flame 
(Niemeyer & Woosley 1997; Khokhlov et al. 1997). Note 
that whether the DDT occurs by this mechanism is contro- 
versial (Niemeyer 1999), and thus the exact density of the 
DDT is stfll debated (Niemeyer & Kerstein 1997). There- 
fore, nucleosynthesis constraints on the DDT density are 
important to obtain. 

Motivated by the discussion above we transform the slow 
deflagrations WS15 and CS15 artiflcially into detonations 
when the density ahead of the flame decreases to 3.0, 2.2, 
and 1.7 x 10^ g cm~3 (DD3, DD2, and DDI, respectively, 



where 3, 2, and 1 indicate p7=p/10^ g cm^^ at the DDT). 
Then the carbon detonation propagates through the layers 
with p < 10® g cm"'^. Figures g and ^ show the density 
distribution and expansion velocity after the passage of 
the slow deflagration and the subsequent delayed detona- 
tion. The explosion energy E and the mass of synthesized 
s^Ni of these WSDD/CSDD models as well as W7 and 
W70 are summarized in Table 1. Here W70 is the same 
hydrodynamical model as W7 except for the initial mass 
fractions of ^(^^Ne) = 0.0, X{^^C) = 0.50, Xi^'^O) = 
0.50. This corresponds to zero initial metallicity because 
^^Ne originates from the initial CNO elements. 

The increase of the ^^Ni mass from W7 to W70 is due 
to the fact that the composition in W70 corresponds to 
a ye=0.5 or a proton/neutron ratio of 1, i.e., symmetric 
matter. In mass zones that are not affected by electron 
capture but which undergo complete Si burning, ^^Ni is 
then the dominant nucleus, without competition by more 
neutron- rich species. Detonations in lower density matter 
do not necessarily lead to complete Si burning. There- 
fore, the amount of ^^Ni is smallest in DDI models and 
largest in DD3 models. The small difference between WS 
and CS DD-models is again a Ye effect. CS models have a 
smaller central density, which leads in the central regions 
to a smaller number of electron captures and larger 
values, resulting again in more symmetric matter. 

2.4. Input Physics 

The methods for performing these calculations were the 
same as used in Shigeyama et al. (1992) and Yamaoka 
et al. (1992). We apply an implicit Lagrangian hydro- 
dynamics code (Nomoto et al. 1984) for the slow defla- 
gration and a Lagrangian PPM hydro code, as used in 
Shigeyama et al. (1992), for the detonation phase. Both 
codes use the same mass grid of 200 radial zones for the 
white dwarf. The nuclear reaction network and the re- 
action rate library utilized are the same as described in 
Thielemann et al. (1996), i.e., thermonuclear rates using 
the Hauser-Feshbach formalism (Thielemann, Arnould, & 
Truran 1987), experimental charged particle rates from 
Caughlan & Fowler (1988), neutron induced rates from 
Bao & Kappeler (1987), and extensions towards the pro- 
ton and neutron drip lines from van Wormer et al. (1994) 
and Rauscher et al. (1994) with Coulomb enhancement 
factors from Ichimaru (1993). Electron capture rates were 
adopted from Fuller, Fowler, & Newman (1980, 1982, 
1985). For nuclei beyond ^=60, only ground state de- 
cay properties were used. We consulted the analysis of 
Aufderheide et al. (1994) so that the influence of a nu- 
cleus with a signiflcant impact on Yf, (either via decay or 
electron capture) was neglected in either of the conditions 
experienced in our calculations. It might, however, require 
a further study to test against the sensitivity of these elec- 
tron capture rates. The present set of Fuller, Fowler & 
Newman (1982) is based on estimates for average proper- 
ties of the Gamow- Teller giant resonance rather than on 
more secure shell model calculations for fp-shell nuclei in 
the Fe-group (Dean et al. 1998). 

One of the problems of nucleosynthesis calculations 
which follow a thermonuclear evolution on long timescales 
through high-temperature regimes is the lack of accu- 
racy. This lack of accuracy is due to the fact that in 
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such situations, where actuaUy a nuclear statistical equi- 
librium should exist, the cancellation of huge opposing 
rate flows is only attained up to machine accuracy (i.e., 
« 10~^^ — 10~^^x the term size). In a similar way the 
term 1/Af, appearing in the diagonal of the Jacobi matrix 
of the multidimensional Newton- Raphson iteration (Hix & 
Thielemann 1999b), can become numerically negligible in 
comparison to reaction rate terms in the same sum, which 
leads to numerically singular matrices. For that reason, 
we have adopted here the accurate solution, i.e., we fol- 
lowed the nuclear evolution with a screened NSE network 
containing 299 species during periods when temperatures 
beyond T=6 x 10^ K are attained (see Hix & Thielemann 
1996, 1999a). This takes into account the changes in bind- 
ing energy or reaction Q-values due to screening. Weak 
reaction rates (electron captures and beta-decays), which 
are not in an equilibrium and occur on longer timescales, 
were included by using the NSE abundances. In this way 
we could track accurately the evolution of the electron 
fraction Ye- 

3. EXPLOSIVE NUCLEOSYNTHESIS 
3.1. Slow Deflagration 

In the present study we follow the approach discussed 
above by varying the ignition density and the initial defla- 
gration velocity and test specifically the effect on the Fe- 
group composition in the central part. Figure ^ shows the 
maximum densities and temperatures for the inner mass 
zones experiencing complete and incomplete Si burning in 
the models WS15, WS30, CS15, and W7 as a comparison. 
It can be recognized that WS15 and WS30 experience den- 
sities similar to that of W7, which is due to the common 
ignition densities of initial model W, with WS15 attaining 
the highest temperatures. CS15, based on initial model C, 
is shifted to smaller densities. 

During the burning the central region undergoes elec- 
tron captures on free protons and iron peak nuclei. The 
central densities, similar to W7 for both WS models but 
combined with higher temperatures, result in more en- 
ergetic Fermi distributions of electrons and larger abun- 
dances of free protons. This leads to larger amounts of 
electron captures on free protons and nuclei and smaller 
central Ye values. Figure || shows the final Ye during 
charged-particle freeze-out and before long-term decay of 
unstable burning products. Table 2 lists this value of 
in the center, Ycc, together with nucleosynthesis informa- 
tion, which will be discussed later. In summary, y^c is 
lower for higher central densities and slower deflagrations, 
the latter leading also to higher central temperatures. Fig- 
ure I and Table 2 also list two models WSL and WLAM 
that experience (partially or entirely) deflagration veloc- 
ities as small as the laminar speed, which apparently do 
not follow this logic anymore. All models are discussed in 
detail in the following paragraphs. 

WS15, corresponding to a slow deflagration with a 
burning-front propagation of 1.5% of the sound speed, 
reaches higher central temperatures than WS30. The case 
WS30, corresponding to a burning- front propagation of 
3% of the sound speed, has similar densities, but slightly 
lower temperatures for each radial mass zone in compar- 
ison to WS15 (also for a shorter duration). Besides the 
different central value, both models also have a different 



central Y^ gradient (see Fig. ||). The WS15 curve, with 
lower central values, reaches ye=0.4985 (inherited from 
the progenitor white dwarf after He burning) at smaller 
radii than WS30. A smaller deflagration speed causes a 
later arrival of the burning front at a given mass coordi- 
nate. Thus, matter at this mass coordinate has a longer 
time to preexpand between the arrival of the central in- 
formation with sound speed and the arrival of the burning 
front. Therefore, burning occurs there at a lower den- 
sity (and temperature) with smaller average electron en- 
ergies, causing less electron capture. This means that a 
smaller propagation speed produces a smaller central Ye, 
but reaches Ye — 0.4985 also at smaller radii, i.e., produces 
a steeper Ye gradient. W7, with a different description of 
the burning-front velocity, but on average a larger speed, 
led to a higher central Ye and a flatter Ye gradient. 

CS15 and CS30, which also corresponds to a burning 
front propagation with 1.5% or 3% of sound speed, be- 
have in a fashion very similar to WS15 and WS30 but are 
characterized by a smaller central ignition density. Be- 
cause of a smaller ignition density the central Ye values are 
larger. Otherwise the Ye-gradients are the same for CS15 
and WS15 as well as CS30 and WS30, each pair having 
the same burning-front speed. CS50, a case with 5% of 
sound speed, has a more extended region of decreased Ye 
out to larger masses (coming close to the behavior of W7), 
but because of the lower ignition density and lower central 
temperatures the central Y^^c is larger again. 

WSL is an additional case and corresponds to a calcula- 
tion with an initially laminar flame speed (out to 0.05 Mq, 
see Fig. |]), before an artificial acceleration is induced 
via turbulent mixing of ashes. WSLAM starts in exactly 
the same manner; however, the front stays laminar, i.e., 
only the minimum fiame speed is permitted. According to 
the previous discussion one would expect an even steeper 
Yg-gradient, and thus a lower central Yg, owing to higher 
temperatures during a longer duration time. The steeper 
gradient can be seen outside 0.06 Mq for WSL, where the 
front starts to accelerate beyond the laminar speed. Inside 
the range of 0.06 MQthe Ye is almost constant and larger 
than e.g., in WS15, which experienced the same ignition 
density and larger flame speeds. There are two reasons for 
this behavior in the very central zones: 

1. Initially during the burning, the lowest central Yg- 
values of ~ 0.432 are obtained as expected, but this value 
is - in contrast to all other cases encountered previously 
- smaller (more neutron-rich) than the Fe-group nuclei in 
the valley of stability (see the most neutron-rich entries in 
Table 2, i.e., 5°Ti, ^^Cr, ^^Fe, and ^^Ni with Z/A values 
of 0.44, 0.444, 0.448, and 0.438). This leads to counterbal- 
ancing /3~ -decays which win against electron captures as 
densities and temperatures (slowly) decrease while matter 
is still in an NSE. This causes an increase in Ye to the dis- 
played values before charged-particle freeze-out (see Fig. 

2. During the laminar front propagation, the large 
amount of central electron captures leads to neutrino losses 
which reduce the local energy release to a point where the 
expansion is very small (see the time dependence of p and 
T in Fig. ^). This keeps temperatures and densities high 
for a prolonged period, and matter stays in an NSE (the 
differences between WSL and WLAM emerge when the 
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burning front is accelerated in WSL beyond the laminar 
speed). This NSE distribution of nuclei - with a total 
more neutron-rich than stable Fe-group nuclei - permits 
/3~ -decay of short-lived nuclei toward a total corre- 
sponding to neutron-rich, stable Fe-group muclei (beyond 
0.44) as long as temperatures are high enough to ensure 
an NSE. This results in the surprising fact that a burning 
front that stays laminar (WLAM) causes a higher Yg af- 
ter charged-particle freeze-out than WSL (see changes in 
Fig. H after 20 s). Thus, a very small Ye, which produces 
e.g., large amounts of '^^Ca (Z/A=0.417), can only be at- 
tained in models with a higher ignition density, causing 
large amounts of electron capture, but with nonlaminar 
burning fronts that permit a fast expansion rather than 
keeping matter in NSE for a long time. 

With the exception of the laminar burning-front models 
ust discussed, which show a different behavior in Figure 
J, we can summarize the major options for slow deflagra- 
tions to change Yg values in the central part of SNe la 
explosions: (1) the burning- front speed determines the Yg 
gradient, and the slowest speeds lead to the smallest cen- 
tral values; (2) lower central ignition densities cause larger 
Ye values, with the gradient, however, depending only on 
the propagation speed. These features hold true as long 
as Ye values in explosive burning do not drop below 0.44, 
when competing /3~ -decays have also to be taken into ac- 
count. While the correct conditions occurring in SNe la 
might be still forthcoming from detailed multidimensional 
hydrodynamic calculations, this parameter study shows 
how nucleosynthesis results can give important clues to 
Vdcf and Pc^ign- 

Figures ||ab, ^b, [To|ab, and p^b show abundance plots 
(mostly of Fe-group nuclei) for the central parts of CS15, 
WS15, CS30, WS30, CS50, W7, WSL, and WLAM where 
electron capture plays an important role. We see that 
Ye values of 0.47-0.485 lead to dominant abundances of 
^''Fe and ^*Ni; values between 0.46 and 0.47 produce dom- 
inantly ^^Fe; values in the range of 0.45 and below are 
responsible for ^^Fe, ^^Cr, 5°Ti, and ^'^^^'^Ni; and values 
below 0.43-0.42 are responsible for ''^Ca. Figure || clearly 
indicates that because of the flat Ye-gradient of W7, the 
total amount of matter experiencing the range of Ye=0.47- 
0.485 is much larger than in cases with slower deflagration 
speeds and larger Yg-gradients. WS15 and CS15 have sim- 
ilar total amounts of ^''Fe and ^*Ni, but at slightly differ- 
ent locations (see the Yg-range in Fig. || and the different 
mass scales in Figs. ^ and ||b). WS30 and CS30 con- 
tain a larger amount of ^''Fe and '^^Ni, owing to the flatter 
Ye-gradient. CS50 comes close to the original W7 model. 

3.2. Fe-Group Composition in Slow Deflagrations 

One of the motivations for the present exercise is to get 
constraints from comparison with solar Fe-group abun- 
dances. This can provide a better understanding of the 
burning-front propagation and test how the otherwise 
quite compelling features of the widely used model W7 
can be improved. Figure shows the ratio of abundances 
produced in W7 to solar abundances. These are the re- 
sults of recalculations of the original W7 with the present 
reaction rate library and an increased accuracy in mass 
conservation in comparison to earlier studies due to the 
screened NSE treatment at high temperatures discussed 



in §2.4. Displayed are abundance ratios after the decay 
of unstable nuclei, normalized to unity for ^^Fe. If SN la 
events are a relatively homogeneous class, the comparison 
of nucleosynthesis products with solar abundances is actu- 
ally meaningful without averaging over a complete sample. 

It is immediately obvious from Figure ^ that the pro- 
duction of Fe-group nuclei in comparison to their solar 
values is a factor of 2-3 larger than the production of in- 
termediate nuclei from Si to Ca. When considering that 
SNe la produce about 0.8Mq of Fe-group nuclei in com- 
parison to wO.IMq from SNe II, and that the Ia/(II-|-Ibc) 
ratio is about 0.15-0.27 in our galaxy (van den Bergh & 
Tammann 1991; Cappelaro et al. 1997), SNe la are respon- 
sible for more than 55% of Fe-group nuclei. Thus, even if 
an isotope has no contribution from SNe II, this implies 
that the isotopic ratios among the Fe-group in the SNe la 
ejecta should not exceed the solar ratios by a factor of ~ 2, 
in order to result in solar ratios for SNe II -t- SNe la. If 
we assume on the other hand that an (Fe-group) isotope 
is made by SNe II in solar proportions, then it has also to 
be made in solar proportions in SNe la to obtain the solar 
mix for the sum of both contributions. An overproduction 
by SNe II would even ask for an underproduction in SNe 
la. If we are conservative and neglect the latter case, an 
overproduction of a factor 1-2 in SNe la would be permit- 
ted. Multiplying this with an average uncertainty factor 
of 2 would permit overproductions of 2-4, dependent on 
the fact whether an SNe II contribution is existing or not. 
In general we do not know this at present, and we take an 
overproduction factor in SNe la of ~ 3 as an alarm sign for 
nucleosynthesis constraints. In this respect, we notice in 
Figure 12 the overproduction of ^"^Cr and ^*Ni by a factor 
of ^ 4 — 5. Here ^''Cr is an = Z + 6 nucleus originat- 
ing from the very central regions with low Ye, while '^^Ni 
(and ^^Fe) are nuclei with N = Z + 2 measuring the bulk 
neutron excess of the material affected by the deflagration 
wave in the intermediate Ye range ^ 0.48 (see also Table 
2). Outside of 0.3Mq in the exploding white dwarf, where 
electron capture is not effective, the neutron excess is only 
determined by the ^^Ne {N = Z + 2) admixture to ^^C and 
^^O in the original composition, stemming from ^'*N in He 
burning, which in turn originated from all CNO-nuclei in H 
burning. The Yg or neutron-excess t] outside O.3M0 is thus 
a measure of the metal abundance (nuclei heavier than He) 
and the galactic age of the white dwarf. The quoted cal- 
culations were performed with Ar(^^Ne) — 0.025, which 
corresponds to of ~ 30 % higher than the solar metallic- 
ity, thus overestimating the average value for "old" white 
dwarfs undergoing a SN la event. 

Using an averaged metallicity of SN la progenitor sys- 
tems less than solar would reduce the overproduction of 
these nuclei. This is shown in Figures 13| for a deflagra- 
tion model like W7 but with zero metallicity (W70). Fig- 



14 1 and 14b give the comparison of the abundance 



distributions. Bravo et al. (1992) have performed a sim- 
ilar test. As the total SNe la contribution in our galaxy 
is given by an integral over time or metallicity up to the 
formation of the solar system, one expects average metal- 
hcities of 0.5-0.6 times solar. This brings the ^'^Fe/^^Fe, 
57pg/56pg^ ssNi/s^Fe, and ^^Ni/sepe ratios within a fac- 
tor of 2-3 of solar, respectively, which corresponds to the 
present uncertainty range of thermonuclear reaction rates 



7 



and the minimum sa50% contribution of SNe la to the Fe 
group. It does not reduce the ^''Cr abundance sufficiently, 
which originates solely from the central layers, where Ye is 
the smallest and entirely due to electron captures rather 
than the metallicity in terms of ^^Ne. An interesting as- 
pect of the change in metallicity, leading to reductions of 
^''Fe and ^*Ni, is the varying amount of early Fe (^''Fe 
before ^^Ni decay) and late Ni (^*Ni after ^^Ni decay) in 
SN la ejecta, leading to features that can be analyzed in 
observed spectra (see Hoflich et al. 1998). This analysis 
covers one uncertainty regarding the initial composition. 
Another one would be a variation of the initial 
ratio, depending strongly on the initial white dwarf mass 
and metallicity. This has also been addressed recently by 
Umeda et al. (1999) and Hoflich et al. (1999). Thus, while 
the metallicity and initial composition are one set of pa- 
rameters, the ignition density and burning-front velocity 
represent another set, as outlined in §2. A burning front 
with a smaller velocity could reduce the amount of mate- 
rial in the range 0.47-0.485, where '^^Fe, ^^^Ni, and ^"^Ni 
are produced in large amounts (see e.g.. Table 2). 

In Figures [Tsl a and 15b, ^6|a and 16b, |l^, and ^8^ and 
18b the ratios to solar abundances (normalized to ^^Fe) 
are displayed. Here the results of the central slow defla- 
gration studies have been merged with (fast deflagration) 
W7 compositions for the outer layers, where Y^ is given 
by the initial ^^Ne and not by electron captures. Thus, 
these are not yet full delayed detonation models (which 
will be discussed in the next subsection), but more pre- 
liminary approximations to test the central Fe-group re- 
sults. In comparison to Figure ^ we see that in all cases 
the ^*Ni problem is strongly reduced and would be fully 
resolved when using also smaller metallicities (see W7 vs. 
W70). This is due to the steeper Y^ gradient which pro- 
duces less matter in the intermediate Ye range 0.47-0.485. 
The smaller propagation speed has, however, also the con- 
sequence that Ye dips deeper in the central layers of WS15 
and WS30 than in W7. Such Ye's lead to the overproduc- 
tion of 5°Ti and ^^Cr. This is not the case for CS15 and 
CS30, due to the lower ignition density. 

Another interesting point surfaces, which was also ad- 
dressed preliminarily in Thielemann et al. (1996) and also 
Meyer, Krishnan, & Clayton (1996) and Woosley (1997b). 
If one takes the results of presently existing and still crude 
SNe II nucleosynthesis calculations from initiated explo- 
sions as well as the results of W7, it turns out that for 
some intermediate mass and all Fe-group elements the 
most neutron-rich nuclei are drastically underproduced. 
The central Ye values of slow deflagration models comes 
close to conditions, where these nuclei are produced in a 
normal freeze-out. This has the effect that nuclei like ^°Ti, 
54Cr, sspe and partially ^^Ni or ^SCa are produced for Ye 
values below 0.46 (or 77 = 1 — 2Ye = 0.1). Whether this 
leads after integration over all mass zones just to solar 
abundances or to a strong overproduction will be tested 
here (see also Table 2). Khokhlov et al. (1992) under- 
took already a preliminary assessment of this question but 
did not consider all of these nuclei in their calculations. 
Woosley (1997b) tested models with higher ignition den- 
sities and burning-front velocities, whether such elements 
can be produced. Here we take again the philosophy to 
obtain constraints for the " average" case of SNe la, to test 



on the one hand whether they can produce such isotopes 
at all, and if so, whether constraints can be set for the 
conditions in order to avoid overproductions beyond the 
often mentioned factor of 3. 

It is a well-known fact that SNe la are not identical (see 
e.g., Hamuy et al. 1995) and that therefore a continuous 
superposition of such models has to be responsible for the 
solar Fe-group composition. On the other hand, we also 
know that the majority of SNe la come from a narrow 
window of conditions (Branch 1998), i.e., the notion of an 
average SN la event makes sense. If Fe-group elements 
are produced by SNe la to ~ 55 % and la's were the only 
sources of the neutron-rich Fe-group nuclei, we must con- 
clude that an overproduction of a factor of 3 is permitted 
and that a larger overproduction has to be avoided for the 
average event. We see that for the cases discussed before 
(CS15, WS15, CS30, WS30) ^^Fe and ^^Ni can be pro- 
duced within the permitted uncertainty limits C^^Ni how- 
ever having a tendency for overpoduction). These are the 
nuclei with intermediate Ye (0.47-0.485), which measure 
the Ye gradient via the amount of mass contained in that 
interval. CS50, which has a larger burning-front speed and 
a flatter slope, starts to overproduce these nuclei. 

The more neutron-rich Fe-group nuclei depend more on 
the central layers, which experience lower Yg-values. ^^Cr 
is produced within permitted limits in CS15 and CS30 but 
is overproduced in WS15 and WS30 owing to the higher 
ignition densities and central Ye values that are too low. 
^"Ti is produced close to solar values for CS15, under- 
produced by CS30, well produced for WS30, and clearly 
overproduced in WS15, preferring the central Ye values 
of CS15 and WS30, which are similar. ^''Ni is essentially 
only made in large quantities for such low- Ye conditions 
as in WS15. Close to solar values cannot be attained for 
*®Ca in any of these models, and its production would re- 
quire lower Ye values. An attempt to do this was a purely 
laminar, i.e., the slowest possible, burning front. The dis- 
play in Figure |l8| a and 18b makes clear that this is not 
possible, as discussed in §3.1. The main reason is that 
such conditions produce nuclei which are unstable against 
P~-decay. In intermediate phases sufficiently low Ye val- 
ues are attained as a result of electron captures. NSE 
or Quasi-statistical equilibrium (QSE) redistributes abun- 
dances according to the then obtained smaller Ye values. 
/3~-decay of the short-lived neutron-rich isotopes leads to 
an increase in Ye during the expansion, when the densi- 
ties and temperatures (and therefore the electron capture 
rates) decrease, before freeze-out from charged particle re- 
actions and NSE/QSE (see Fig. ^). Thus, to produce a 
nucleus like *^Ca in sufficient amounts, only very high- 
density ignitions (followed by a fast expansion to avoid 
the influence of /3~-decays) of progenitor systems, which 
barely avoid an accretion induced collapse (AIC), might 
be responsible (see, e.g., Nomoto & Kondo 1991; Woosley 
1997b). 

From this exercise we see that not all of these neutron- 
rich nuclei can be made in similar proportions for one set 
of deflagration parameters unless a detailed fine-tuning of 
^dcfif) and Ye(r) is performed or multidimensional prop- 
agation of the burning front produces exactly a superpo- 
sition of conditions as needed to fit all abundance con- 
straints. In general CS15 and CS30 seem to be better 
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models than WS15 and WS30 in terms of avoiding a large 
overproduction of neutron-rich elements. The burning- 
front speed in the central layers seems to be constrained 
to values below 5% of the sound speed in order to avoid 
the overproduction of ^^Ni (see CS50 and Figure ^ as 
well as Table 2 with comparable problems found for W7 
and WTO). If on the other hand there were no other sites 
than SNe la to produce isotopes like ^°Ti, ^^Cr, ^Spe, 64Ni 
or ^^Ca, we would have to overproduce in comparison to 
solar by about a factor of 2, and some features of the 
models WS15, WS30 and even higher density events were 
needed, which can fill in the remaining deficiencies. Thus, 
one would need the majority of events similar to CS15 or 
CS30 and a smaller number of higher density ignitions to 
produce the more neutron-rich nuclei. This could guar- 
antee on the one hand some overproduction, which (in 
combination with SNe II) would permit solar abundances 
in total, and avoid on the other hand not-permitted over- 
production. 

3.3. Delayed Detonation 

The final aim, after constraining the central slow defla- 
gration part of SNe la, is to find also composition con- 
straints for the deflagration detonation transition (DDT). 
As discussed in §2.3, we chose transition densities of 3.0, 
2.2, and 1.7x10'' g cm~^ for the given models WS15 
and CS15, which turn them into WS15DD3, WSlSDm 
WS15DD1 or CS15DD2 and CS15DD1. Figures |l^ - ^ 
show the abundance distributions of slow deflagrations 
combined with delayed detonation models against the ex- 
pansion velocity and Mr of DD models. The central re- 
gions of these models have been shown before in Figures 
H and ||, thus the abundance distributions of neutron-rich 
species such as ^''Cr, ^°Ti, ^^Fe, and ^^Ni are not repeated 
here. 

As the deflagration wave propagates outward, the white 
dwarf gradually expands to undergo less electron capture 
and thus mostly ^^Ni is synthesized. Eventually, the de- 
flagration enters the region of incomplete Si burning and 
explosive 0-Ne-C-burning, where the transition to a det- 
onation occurs. For comparison the abundance distribu- 
tions of W7 and WTO were shown in Fugure Ol The total 
masses of ^^Ni(^^Fe) produced in these combined models 
have been summarized in Table 1. 

These theoretical abundance distributions can be com- 
pared with the observed expansion velocities of several el- 
ements as estimated from supernova spectra. It is seen 
that WDD2 and WDDl produce two Si-S-Ar peaks at low 
velocity (~ 4,000 km s~^) and high velocities (10,000 - 
15,000 km s~^). The intermediate mass elements at low 
velocities are important to observe at late times in order to 
distinguish between models. In particular, the minimum 
velocity of Ca in WDD models is ~ 4,000 km s~^, which 
would be higher for a faster deflagration. 

The Ca velocities should be compared with the observed 
minimum velocities of Ca indicated by the red edge of the 
Ca II H and K absorption blend (Fisher et al. 1995). The 
lowest velocities of O and Mg also provide interesting con- 
straints. For example, SNe 1990N, 1992A, and 1991T show 
O in the wide velocity range from « 10,000 to 20,000 km 
(Leibundgut et al. 1991a; Jeffery et al. 1992; Mazzali 
et al. 1993; Kirshner et al. 1993). For W7 and WDDs, 



the minimum O velocity is 12,000 - 15,000 km s"^ The 
observed O velocity as low as 10,000 km s^^ may indicate 
a mixing of O in the velocity space. 

Meikle et al. (1996) have observed a P Cyg-like feature 
at - 1.05/1.08 fim in SN 1994D and 1991T. They note 
that, if this feature is due to He, He in SN 1994D is likely 
to be formed in an alpha-rich freeze-out and mixed out 
to the high-velocity layers (~ 12,000 km s~^). The max- 
imum velocity of He is 5,000 - 6,000 km s^^ in WDDs, 
being slower than ~ 9,000 km s~^ in W7, so that more 
extensive mixing of He would be required for WDDs than 
in W7. (Note that ^He is plotted neither in Fig. |lj nor 
in Figs. |l^ - |2^, but its location can be easily identified 
with the region where the ^*Ni abundance forms a plateau 
up to the point where it turns over to ^''Fe on the same 
level. This is the transition from alpha-rich freeze-out to 
incomplete Si burning.) Alternatively, if the feature is due 
to Mg, the Mg velocity is confined to 12,500 - 16,000 km 
s-i in SN 1994D, which is consistent with W7 (13,000 - 
15,000 km s~i). For WDDs, on the other hand, the mini- 
mum velocities of Mg are 14,500 km s'^ (WDDl), 16,500 
km s-i (WDD2), and 18,000 km s~i (WDD3), and the 
latter two models seem to have too high velocities. 

Mixing of ejected material in velocity space could oc- 
cur convectively during the propagation of the deflagra- 
tion wave (Livne 1993). Nonspherical explosions induced 
by delayed detonations could also produce nonspherical 
abundance distribution, i.e., elemental mixing in the ve- 
locity space. From the calculated synthetic spectra and 
their comparison with the observations, more advanced 
methods (Harkness 1991) do not favor mixing opposite to 
initial suggestion by Branch et al. (1985). 

4. YIELDS OF SNE lA AND GALACTIC CHEMICAL 
EVOLUTION 

4.1. Features of SN la Nucleosynthesis 

Complete isotopic compositions of WDD and CDD mod- 
els are given in Table 3. Table 3 assumes full decay of all 
unstable species. We provide separately in Table 4 abun- 
dances of long-lived radioactive nuclei, of importance ei- 
ther for gamma-ray detection, extinct radioactivities, or 
chemical evolution. The abundances are compared with 
solar abundances in Figures ^ - 12^, which are normalized 
to ^^Fe. These Figures complement the earlier Figures [ij 
- |l^, where the fast deflagration W7 was utilized in the 
outer layers rather than delayed detonation models. The 
major conclusions on the Fe-group composition remain the 
same as discussed before, in §3.1, owing to the fact that 
they are given by the central slow deflagrations. What 
can be studied here is the additional variation in the ratio 
of Fe-group to intermediate mass nuclei, which depends 
on the deflagration-detonation transition. Of course the 
variation in ^^Fe (originating from ^^Ni) in the outer det- 
onation layers also influences the ratio of neutron-rich Fe- 
group isotopes (from central locations) to ^^Fe. Table 3 
includes besides all DD-models (CDDl short for CS15DD1 
etc.) also W7 and W70 updated with the latest reaction 
rate set and improved accuracy by using a screened NSE 
treatment for long duration times at temperatures beyond 
6 X 10^ K (Hix & Thielemann 1996). The ratios to solar 
abundances for W7 and W70 were shown in Figures |l^ and 
13. The essential features due to the DDT, as displayed 
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in Figures ^ - |2j and in Table 3 can be summarized as 
follows: 

1 . The synthesized amount of Fe and thus the ratio be- 
tween the intermediate mass elements and Fe, Si-Ca/Fe, 
is sensitive to the transition density from deflagration to 
detonation, as was already shown in Table 1. Among the 
WDD models, WS15DD2 produces only about 25% more 
56Ni than W7 (- 0.6 Mq) but more Si-Ca than W7 by 
40% (Fig. [23|b), since more oxygen is burned in the outer 
layers. Therefore, the Si-Ca/Fe ratios are moved up to a 
certain extent. WS15DD1 has even larger Si-Ca/Fe ratios, 
which are close to solar ratios (Fig. |2^a). This is not indi- 
cated for SNe la, owing to observations of low-metallicity 
stars reflecting the average SN II behavior (a reasoning 
outlined below in more detail). However, direct observa- 
tions of the Si-Ca/Fe ratio in SNe la remnants (Tycho, 
SN 1006, etc.) are highly important and needed in order 
to distinguish between the models (Hughes et al. 1995; 
Miyata et al. 1998; Hwang, Hughes & Petre 1998). 

2. Neutron-rich species such as ^"'Cr and ^°Ti are mostly 
produced in the slow deflagration phase. The degree of 
their overproduction with respect to ^^Fe depends also on 
the mass of ^^Ni produced in the outer detonation layers, 
as seen in Figure ^ This also explains why in Table 2 
the entries for the same central models (WS15, CS15 etc.) 
change for neutron-rich species, although the central part 
from which these neutron-rich species originate is unaf- 
fected by the detonation. The reason is that the ratios in 
comparison to ^^Fe are taken. 

3. There are some Fe-group contributions from alpha- 
rich freeze-out and incomplete Si-burning layers that de- 
pend on the DDT. Figures |l^ - |2^ show that the mass re- 
gion experiencing incomplete Si burning (indicated by the 
^"^Fe plataeu) decreases in the sequence DD1-DD3. The 
region indicated by the ^*Ni plateau experiences alpha- 
rich freeze-out (the He abundance is not shown here) and 
increases in the sequence DD1-DD3. '^^Fe (decaying to the 
dominant Cr isotope ^^Cr) and ^^Co (decaying to the only 
stable Mn isotope ^^Mn) are typical features of incomplete 
Si burning. ^^Cu (decaying to the only stable Co isotope 
^^Co) is a typical feature of an alpha-rich freeze-out. For 
these reasons we see the strongest ^^Cr and ^^Mn over- 
abundances in DDI and the strongest appearance of ^^Co 
(while still underabundant) in DD3. 

4.2. The Role of SN la and SN II Contributions 

The chemical evolution of galaxies is dominated by its 
main contributors SNe II, SNe la, and planetary nebulae. 
The latter do not contribute to the element abundances 
in the range O through Ni (although a few specific minor 
isotopes can be produced in the s-process). Thus for the 
aspects considered here, we have to explain galactic evo- 
lution and also solar abundances by the combined action 
of SNe II and la. The ratio of Fe-group elements to Si-Ca 
in SNe la is of specific importance, in order to see how 
the overabundance of 0-Ca/Fe (in comparison to solar) 
in SNe II can be compensated. Combined nucleosynthe- 
sis products of SNe la and SNe II with varying ratios can 
be compared to solar abundances. An important aspect 
for such an undertaking is, however, to test the individ- 
ual components against existing observations first, before 
trying to attain a good solar mix (with possibly wrong pre- 



dictions for the individual SN I and SN II components). 

Although we do not have a good quantitative measure 
from SN la observations for the Fe-group to Si-Ca ratio, 
the fast deflagration model W7 seems to give a good over- 
all agreement via synthetic spectra calculations with ob- 
served la spectra (see Branch 1998 for a review). W7 has 
speciflc deficiencies in the global isotopic Fe-group com- 
position from the inner layers, as discussed before, e.g., 
with respect to ^^Ni and ^"^Cr, but this does not affect the 
element ratios too strongly. 

A similar or even worse situation is found for the quanti- 
tative analysis of SN II spectra. However, an independent 
tool exists that measures the integrated SN II yields: ob- 
served surface abundances of low-metallicity stars, which 
witness the abundances of the interstellar medium at their 
point of formation during early galactic evolution, when 
only SNe II contributed. For metallicities in the range 
—2 <[Fe/H]< —1 we expect the integrated (mass aver- 
aged) properties of SNe II (see e.g., Nakamura et al. 1999). 
Here [x/y] is defined as logio[(x/y)/(x/y)Q]. Such abun- 
dance features were reproduced with nucleosynthesis prod- 
ucts of SNe II as a function of stellar mass, taken from the 
calculations by Nomoto & Hashimoto (1988), Hashimoto 
et al. (1996), and Thielemann et al. (1996) as summa- 
rized in Hashimoto (1995), Tsujimoto et al. (1995), and 
Nomoto et al. (1997a). SNe II yields, integrated from 
mi = 10 M0 to niu = 50 M© with a Salpeter IMF, 
are also given in Table 3. The upper mass bound m„ is 
chosen to reproduce [0/Fe] = -1-0.4, which is consistent 
with the observations of low metallicity stars for [Fe/H] < 
- 1. Such observations in low-metallicity stars give also 
the best constraints on average Fe-group abundances of 
SNe II, which are poorly known theoretically owing to the 
still existing lack of self-consistent core collapse supernova 
models. The representation of the Fe-group (beyond Ti) 
closest to the low metallicity observations seems to be the 
one of our 2OM0 star (see Figure 5b in Thielemann et al. 
1996 and note that the observed Co abundance has come 
down to about -0.1, i.e., it shows a better agreement with 
the dashed line). 

A possible deviation from solar ratios has to be made 
up by an opposite behavior of SNe la setting in at about 
[Fe/H] = — 1 in order to attain solar values at [Fe/H]=0. 
Fe-group elements for which information is available are 
Ti, Sc, Cr, Mn, Co, Ni (Magain 1987, 1989; Gratton & 
Sneden 1988, 1991; Gratton 1989; Zhao & Magain 1990; 
Nissen et al. 1994). They lead to typical uncertainties 
of 0.1 dex and one finds average SN II values of 0.25, 0, 
-0.1, -0.3, -0.1, -0.1. Taken the typical uncertainty of 0.1 
dex, this leaves Ti and Mn as elements with clear signa- 
tures for a SN II behavior different from solar, which asks 
for the opposite SN la behavior. Cu and Zn start having 
strong s-process contributions. We avoid their discussion 
because of these complications and because they are not a 
clear indication for the required SN la signature. 

Ti is dominated by the isotope "^^Ti (see Table 3), i.e., 
we have to relate the element ratio Ti/Fe to ^^Ti/^^Fe. We 
see that for all models that underproduce Si-Ca (as needed 
for SNe la), i.e., the DD2, DD3 and W7 models, ''^Ti is 
also underproduced by similar amounts. This agrees with 
the observational trend that Ti is dominantly produced by 
SNe II and can be understood from the ^^Cr abundances 
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(decaying to the main Ti isoptope ''*Ti). ^^Ti is only pro- 
duced in a strong alpha-rich freeze-out as it occurs in SNe 
II. The small region of a weak alpha-rich freeze-out, indi- 
cated in Figure U, is not sufficient. 

The tendency is not so clear for Mn. The only stable Mn 
isotope is ^^Mn, typically produced from unstable ^^Co. 
This isotope results (1) from incomplete Si burning with a 
relatively high Yg (e.g., 0.4985; see its production in SNe 
II in Fig. 1 in Thielemann et al. 1996 and the present 
Figs |l^ - ^ in the incomplete Si-burning regions) or (2 
from a somewhat reduced Yg (around 0.49) in complete Si 
burning (see Figs. ^ - 1|). This can be understood within 
the framework of quasi-equilibrium groups in Si burning, 
in our case the Si and the Fe group. Incomplete burning 
leads to a small total abundance in the Fe group in com- 
parison to the Si group. Hix & Thielemann (1996) found 
that in such a case the composition in the Fe group is typi- 
cally more neutron-rich than expected from the global Yg . 
Hence, ^^Co with a Z/A=0A9 is in incomplete burning 
also produced for conditions with a global Ye of 0.4985. 
Both locations (in complete and incomplete burning) can 
be nicely seen in Figure nj, with the additional Yg infor- 
mation taken from FigureW. As discussed before, the same 
is found in the central parts in Figures ^ - ^ and globally 
in Figures 09^" 
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The problem is, however, that with 
the exception of the DDI models and W7, all DD models 
underproduce Mn/Fe in comparison to solar, opposite to 
what is required from average SN II yields inferred from 
low-metallicity observations. It appears that the dominant 
source for Mn is the incomplete Si-burning region, which 
is most extended in the DDI models. 

^^Cr is the dominant nucleus of the Cr isotopes. DD2 
shows a slight overabundance, DDI a stronger overabun- 
dance. ^^Cr originating from ^^Fe is also a nucleus domi- 
nated by incomplete Si burning and therefore this behav- 
ior is understandable. The light underabundance in SNe 
II can thus be compensated by SNe la. Co from SNe II is 
slightly underproduced. All our calculations show an un- 
derabundance in SN la models. It seems that the alpha- 
rich freeze-out nucleus ^^Cu (decaying to the only stable 
Co isotope ^^Co) is never produced sufficiently. However, 
there exists the same problem in SNe II (Nakamura et 
al. 1999) and possibly other (nuclear?) sources might be 
the origin of this behavior. ^^Ni originates from neutron- 
rich central regions and partially from alpha-rich freeze- 
out (where also ^^Ni is produced via ^^Zn decay). We 
have worked hard on our models to avoid an overproduc- 
tion. Thus a slight underproduction in SNe II can easily be 
compensated. One should also have a look at Bravo et al. 
(1993) and Matteucci et al. (1993), who performed chem- 
ical evolution calculations and addressed some of these 
questions with the then existing observational and model 
constraints. 

We finally want to return to neutron-rich Fe-group iso- 
topes and point to a different witness of galactic evolution. 
While typically astronomical observations can give infor- 
mation only about element abundances (with a few excep- 
tions from molecular lines in stars) , there exists one source 
of isotopic information, the so-called isotopic anomalies 
in meteorites. They are usually contained in "inclusions" 
of primitive meteorites consisting of minerals with much 
higher melting temperatures than the surrounding mat- 



ter. This gives some indication that they originate from 
unprocessed " star dust" that survived temperatures in the 
early solar system and can give direct clues about its stel- 
lar origin. This is essentially proven for some SiC grains, 
graphites and diamonds (e.g., Zinner, Tang, & Andrews 
1989; Zinner 1995; Travaglio et al. 1999). 

The nuclei ^®Ca, ^°Ti, ^"^Cr, and ^*Fe, which play a cru- 
cial role in the central layers of SN la explosions, can also 
be found in isotopic anomalies (Lee 1988; Volkening & 
Papanastassiou 1989, 1990; Loss & Lugmair 1990). They 
occur, however, only in Ca-Al-rich inclusions whose history 
is less clear and some chemical processing in the early solar 
nebula has probably happened. Nevertheless, it is inter- 
esting to investigate whether these observed anomalies can 
lead to any indication of their origin. Harper et al. (1990) 
found a correlation with the r-process nucleus ^^Zn. This 
would have pointed toward a SN II origin if the r-process 
originates from SNe II. On the other hand, Ireland (1994) 
found a probe where no correlation with r-process anoma- 
lies existed. Thus, no very clear indication for their SN II 
or SN la origin (or both?) presently exists. In the present 
paper, however, we applied the working hypothesis that 
their source is given by SN la ejecta, which is also indi- 
cated by the Yg-constraints from SN II abundance obser- 
vations (Thielemann ct al. 1996). In case of an additional 
SN II contribution, even more stringent limits would exist. 

4.3. Ratios of SN la to SN II Events in the Galaxy 

We have full results for nucleosynthesis products of SNe 
la from the four models, WS15DD1-3, CS15DD1-2^W7, 
and WTO, and the results are shown in Figures E3- 
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12| - |13|. CS30 and 50 as well as WS30 were calculated 
only for the central slow deflagration layers and not full 
delayed detonation models. This, as it affects the compo- 
sition within the Fe group, might not be a sufficient quan- 
titative basis with respect to relative abundances within 
the Fe group. But as the central deflagration speeds have 
minor influence on the total abundance of ^^Ni (^^Fe), the 
different DD models are probably sufficient to determine 
the best Si-Ca/Fe ratios. The observational data by Graf- 
ton & Sneden (1991) and Nissen et al. (1994) for [x/Fe] 
at low metallicities (—2 <[Fc/H]— < 1), x standing for 
elements from O through Ca, show an enhancement of the 
alpha elements (O through Ca) by a factor of 2-3 (0.3 to 
0.5 dex in [x/Fe]) in comparison to Fe. This is the clear 
fingerprint of the exclusive contribution of SNe II in early 
galactic evolution. It has long been the aim of chemical 
evolution calculations to explain this behavior, among the 
most recent ones being Tsujimoto et al. (1995), Timmes, 
Woosley, & Weaver (1995), Pagel & Tautvaisine (1995, 
1997), and Kobayashi et al. (1998). 

Tsujimoto et al. (1995) tried to determine the ratio of 
SN la to SNe Ib+H events in the Galaxy by aiming for a 
best fit to an overal solar abundance pattern from O to Ni. 
They made use of the mass-averaged SN II yields from 10 
to 50 Mq as shown in Table 3 and earlier results for W7. 
With the aid of a chemical evolution model they obtained 
a ratio Nii^/Nu = 0.12 of the total number of SNe la to 
SNe II (-l-Ib) that occurred in our Galaxy. This resulted 
from an overal best fit to the observed abundances and 
is consistent with the fact that the observed estimate of 
the SNe la frequency is as low as 10 % of the total su- 



11 



pernova rate. The observational constraints for this ratio 
range from values around 0.15 (van don Borgh & Tam- 
mann 1991) to more recent determinations of 0.27 with 
errors of almost a factor of 2 (Cappelaro et al. 1997). 

Rather than performing a full galactic evolution model, 
we want to concentrate on a typical abundance ratio, 
Si/Fe, and its contribution from SNe II and la. The aver- 
aged SNe II yields from our theoretical models summarized 
in Table 3 correspond to a [Si/Fe] =0.347. The results of 
the 20 M0 model of Thielemann ct al. (1996) would corre- 
spond to [Si/Fe]=0.304, which is close to the value of Graf- 
ton & Sneden (1991) and Ryan, Norris, & Beers (1996) of 
0.3. One should, however, be aware of the typical uncer- 
tainties of 0.1 dex. Varying the SN la models among the 
above given hst (W7, W70, DDI, DD2, and DD3) and 
making use of Si and Fe from the averaged SNe II ejecta 
or the 20 Mq star, leads also to a required Ia/Ib-|-II ra- 
tio in order to obtain the solar mix of Si/Fe with both 
contributions. This ratio is derived in the following way: 



M(Si) 
M(Fe; 





iVia 
Nil 



iViaMia(Si) + 7VnM„(Si) 
iViaMia(Fe) + 7ViiM„(Fe), 
Mii(Si) - M(Si)/M(Fe)0Mii(Fe) 
M(Si)/M(Fe)0Mia(Fe) - Mia(Si) 



(3) 



(4) 



The results for the obtained ratios are shown in Table 

5. It is immcdiatcy obvious that the DDI models do not 
result in a physical solution. With their Si-Ca/Fe ratios 
being almost solar, they cannot compensate for the SNe II 
contribution, which have larger than solar values. There 
is a clear need for SNe la to produce ejecta with Si-Ca/Fe 
less than solar. Otherwise all models seem to be covered by 
the error bars for the la/ (IbH-II) ratios laid out by van den 
Bergh & Tammann (1991) and Cappelaro et al. (1997). 
The DD3 values, especially for the 20 Mq Si/Fe ratios 
being closest to the low-metallicity [Si/Fe], might some- 
what indicate the lower limit. It should be noted that 
although the Si/Fe ratios in the DDI - DD3 models show 
a monotonous behavior (i.e., decHne), this does not lead to 
a monotonous variation in Ia/(Ib+II) ratios, owing to the 
form of equation (4), where differences and not the Si/Fe 
ratios enter. This results in negative values for DDI (due 
to a negative denominator) and goes through a maximum 
(when the denominator starts to be positive but is small) 
before decreasing to the DD3 values with an increasing de- 
nominator. This is due to the fact that Si masses decrease 
with the DDI - DD3 sequence, while Fe increases. 

With a low SN la frequency of Ia/(II+Ib)=0.15 and W7 
abundances, "''^Fe from SNe la amounts to about 55% of 
total ■''''Fe (see the discussion in §3.2). Larger values would 
increase this contribution. When assuming that SNe II do 
not produce neutron-rich species like ^^Cr and ^'^Ti at all, 
we permitted for the abundance ratios between such nuclei 
and ^®Fe an overproduction factor of ~ 4. With the larger 
SN la frequencies of Table 5 and Cappelaro et al. (1997) 
such permitted overproductions would need to be further 
reduced and more stringent limits would apply. 

5. SUMMARY AND OUTLOOK 

From the very early days of explosive nucleosynthesis 
calculations, when no direct connection to astrophysical 



sites was possible yet, it was noticed (Trimble 1975) that 
the solar Fc-group composition could be reproduced with 
a superposition of matter from explosive Si burning with 
about 90% originating from a ¥"6=0. 499 source and 10% 
from a ye=0.46 source. We have discussed in some detail 
in the present paper that the central part of SNe la could 
be this second source, while Thielemann et al. (1996) have 
shown that SN II ejecta with Fe<0.498 could cause serious 
problems. 

New results by Hachisu et al. (1999a, 1999b; see also Li 
& van den Heuvel 1997), which include wind losses in the 
interaction of binary systems, come to the conclusion that 
the majority of SNe la progenitor systems experience hy- 
drogen accretion on white dwarfs at a rate that has them 
grow toward the Chandrasekhar mass through steady El- 
and He-burning. This leads to the single degenerate Chan- 
drasekhar mass scenario (SD/Ch). Binary systems with 
steady H-burning on accreting white dwarfs and effective 
accretion rates as high as M > 1 x 10^^ Mq yr^^ might 
correspond to observed supersoft X-ray sources. This also 
would lead to low central ignition densities of the Chan- 
drasekhar mass white dwarf at the thermonuclear runaway 
(< 2 X 10^ g cm~'^), which correspond to the C series of 
the models discussed in the present paper. A small frac- 
tion can deviate from steady H burning and would experi- 
ence weak hydrogen flashes near the end of the accretion 
history. Such cases would correspond to our W series of 
models and even higher ignition densities. Kobayashi et 
al. (1998) discussed metallicity effects and their influence 
on the delay time for the appearance of SNe la in galactic 
evolution. 

Our nucleosynthesis results of the present paper are 
quite consistent with those scenarios and favor case C. The 
reason is that too high ignition densities lead to a high de- 
gree of electron captures and small Ye values, which would 
cause an overproduction of "''''Cr and ^°Ti in excess of what 
is permitted for SNe la in a solar mix. We have to make 
two reservations here. New shell model calculations (Dean 
et al. 1998; Caurier et al. 1999) indicate that the electron 
capture rates could be substantially reduced in compari- 
son to the rates by Fuller et al. (1980,1982,1985) employed 
here. The application of these new capture rates might 
also permit models of our W series without serious devia- 
tions from the allowed central Ye values. In addition, such 
conclusions deal only with the dominant or average SN la 
events. If more ne\itron-rich nuclei like '**^Ca are also the 
result of SN la nucleosynthesis, an occasional event with a 
significantly higher ignition density is required, a fact also 
consistent with the above scenario. 

Our nucleosynthesis results imply that for the domi- 
nant events the central density of the Chandrasekhar mass 
white dwarf at thermonuclear runaway must be as low as 
or lower than ^2 x 10^ g cm~^, though the exact con- 
straint depends somewhat on the flame speed. Here is 
the point where nucleosynthesis predictions can also be 
of help for providing constraints to the supernova mod- 
eling and the burning-front velocity. A carbon deflagra- 
tion wave, propagating as slow as v^ei/vs ~ 0.015 or even 
slightly slower, would be the ideal choice for the neutron- 
rich species such as '''^Cr and ^''Ti; see cases CS15 and 
WS30. The latter case makes also clear that slightly larger 
ignition densities than 1.5 x 10^ g cm~^ (case C) are per- 
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mitted if the flame speed is increased appropriately. An 
acceleration of the flame specid in the outer part of the 
central layers is permitted to about 3% of the sound speed 
or slightly more, but should clearly stay below 5%. For 
such fast burning fronts the Ye-gradient becomes very flat 
and too much material in the range Ye = 0.47-0.485 is 
produced. This leads to dominant abundances of ^'*Fe and 
^®Ni, a feature that was already prominent in W7, and 
causes an excessive overproduction of ^*Ni in galactic evo- 
lution. Our calculations were performed with a constant 
fraction of the sound speed. An acceleration of the burn- 
ing front is expected (Khokhlov 1995; Niemeyer & Woosley 
1997) and future investigations should include such a time 
dependence, which is in agreement with the findings dis- 
cussed above. 

Finally, nucleosynthesis can also give clues about the 
deflagration-detonation transition (DDT). The most obvi- 
ous consequence of choosing different transition densities 
is the amount of ^^Ni produced in a SN la event. Hoflich 
and Khokhlov (1996) flnd from light curve modeling and 
spectra that the typical ^^Ni mass should be in the range 
0.5-0.7 Mq. This agrees with W7. Among the DD models 
it would ask for a value somewhere between DDI and DD2 
(closer to DD2). Here 3, 2, and 1 stand for DDTs when 
densities ahead of the flame decrease to 3.0, 2.2, and 1.7 
X 10^ g cm~^ . DDI is excluded for other reasons. The 
amount of Si-Ca in comparison to Fe is too large in DDI 
models in order to compensate for the well-known over- 
production of Si-Ca in SNc II during galactic evolution. 
Si/Fe ratios in SN la models require specific Ia/(II-|-Ib) 
ratios in order to obtain a solar mix combined with SN 
II contributions (see Table 5). DD2 seems to be closest 
to the present observational constraints for this ratio by 
Cappelaro et al. (1997). 

Small DDT densities favor larger amounts of matter that 
experience incomplete Si burning. Low-metallicity con- 
straints require an overproduction of Mn (and Cr) in SNe 
la. These elements are mostly made as '^''Co and ''^Fe (de- 
caying to Mn and Cr), which are favorably produced in 
incomplete Si burning and would also require a DDT be- 
tween DDI and DD2. (One should, however, realize that 
a fast deflagration could possibly simulate this as well - 
see the Mn overproduction in Figure 12 - and on the other 
hand that these numbers would have to be rescaled or 
reinterpreted in multi-D calculations.) Thus, combining 
all requirements on the DDT from total Ni yields, Si/Fe 
and Ia/(II+Ib) ratios, as well as speciflc elements favored 
in incomplete Si-burning, we would argue for a DDT den- 
sity slightly below 2.2 x 10^ g cm~^ , i.e., results be- 
tween DDI and DD2. One should, however, be careful 
with these constraints based on spherically symmetric ap- 
proximations of the burning front. Full three-dimensional 
calculations could possibly produce the required ratio of 
matter from incomplete Si burning and complete Si burn- 
ing with alpha-rich freeze-out in a different realization. 

More extended calculations that make use of the con- 



clusions presented here, including a time dependence of 
^'dci/^^sound, the bcst choicc for the DDT density, and a 
detailed galactic evolution model, replacing our compar- 
isons of only the global yields, will be the next step to 
undertake. That would also have to include further tests 
of the metallicity of the exploding object, a topic which 
has gained in importance with the cosmological interpre- 
tation of high-redshift SNe la. For this reason we repeat 
here in Figure 25a-25c some of the Figures 19-21 with a 
continuation to smaller mass fractions and purely plotted 
as a function of velocity, in order to magnify the behavior 
of the outer layers. The differences in velocities between 
CDD and WDD models are negligibly small in these plots. 

As mentioned before, the ^^Fe (in the outer layers not 
affected by electron captures but only by the neutron ex- 
cess due to the initial metallicity) ranging in velocities up 
to 15,000 - 19,000 km s'^ in the models DD1-DD3, is a 
strong indicator of metallicity (see Fig. 14). The minima 
in S, Ar, and Ca, if observed in spectra, with positions 
between 16.000 and 21,000 km s^^ in the DD1-DD3 mod- 
els, could give further clues on the deflagration-detonation 
transition. And finally, any unburned intermediate mass 
elements at higher velocities would give a clear indication 
of the metallicity of the accreted matter. Our calculations 
and plots include only initial compositions of ^^C, ^^O, 
and ^^Ne. Solar abundances of Ca, Ar, S, Si, or Fe would 
correspond to mass fractions of 1.2 x 10~*, 1.6 x 10~*, 
7.6 x 10-"*, 1.3 x 10'^, or 2.3 x lO'^. 

Thus, any future investigations of very early time spec- 
tra, leading to abundance observations at high velocities, 
would provide strong constraints on the SNe la mechanism 
and the relation to the metallicity of the individual SNe 
la explosion. 
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Fig. 1. — Temperatures and densities for the deflagration wave (WS15) and delayed detonation (DDI) as a function of 
radial mass coordinate Mr- The labels indicate different instances in time, 1-6 correspond to the deflagration phase, 7-11 
to the propagation of the detonation after the deflagration-detonation transition. 
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Fig. 2. — The distribution of temperatures and densities as a function of time during the passage of the slow deflagration 
(WS15) and the subsequent delayed detonation (DDI). The labeled numbers indicate different mass zones. The prominent 
spike, occurring for labels > 5, shows the action of the detonation and its fast propagation (small time differences between 
different mass zones). 
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Fig. 3. — The distribution of expansion velocities after the passage of the slow deflagration and the subsequent delayed 
detonation. 
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Fig. 4. — IKftiximuiii densities /^max cUiJ leiiiperatures Tg^max diiriiig Si-buriiiiig of the central layers obtained during the 



propagation of the burning front. The model Wf is compared to three slow deflagrailPons which start (for a density of 
2.1xl0®g cm~'^) with burning front velocities of l.§% (slow deflagration WS15) and 3% (WS30) of the sound speed and 
for an ignition density of 1.37xl0^g cm"*^ with a burning front velocity of 1.5% of sound (CS15). The crosses (C-det) 
correspond to carbon detonations of sub-Chandrasekhar mass models. They will not be discussed here, but it is obvious 
that the smaller temperatures and densities lead to negligible amounts of electron captures. 




Fig. 5. — Yg, the total proton to nucleon ratio and thus a measure of electron captures on free protons and nuclei, after 
freeze-out of nuclear reactions as a function of radial mass for different models. In general it can be recognized that 
small burning front velocities lead to steep Fe-gradient which flatten with increasing velocities (see the series of models 
CS15, CS30, and CS50 or WS15, WS30, and W7). Lower central ignition densities shift the curves up (C vs. W), but 
the gradient is the same for the same propagation speed. This only changes when the attained via electron captures 
during explosive burning is smaller than for stable Fe-group nuclei. Then, /3~-decay during the expansion to smaller 
temperatures and densities will reverse this effect (see models WSL and WLAM discussed in more detail in the text). 



19 



.46 



„ .45 



.44 



.43 







5 10 15 20 

Time (seconds) 



CS50 

CS30 

CS15 

W7 

WS30 

WS15 

WLAM 

WSL 



25 



30 



Fig. 6. — Yeit) for the central mass zones in CS, WS, WSL, WSLAM and W7 models. We see that Yg attains the lowest 
values for the highest ignition densities and slowest burning front speeds. This tendency continues for the extreme cases 
of laminar burning front models, but the slow expansion in these models permits fast /3~-decays of short-lived nuclei in 
NSE or QSE equilibria before charged-particle freeze-out. 




Fig. 7. — p{t) and T{t) for the central mass zones in WSL and WSLAM models. The bifurctaion between WSL and 
WLAM marks the point where the burning front was accelerated beyond the laminar speed in WSL, resulting in a faster 
expansion. 



20 





Fig. 8. — Abundance plots (mostly of Fc-group nuclei) for the cases CS15 and WS15. Shown are the locations of the 
nuclei overproduced in W7, ^^Fe and "''^Ni, corresponding to values of 0.47-0.485. Due to the Fe-gradients which are 
steeper than for W7, the amount of matter in a given Fg-range is reduced, but also smaller central values are attained, 
giving rise to more neutron-rich nuclei. A Yg of about 0.46«26/56 (which was also attained in W7) causes no problems 
and leads to a large abundance of stable ^^Fe (not from ^^Ni decay). Values in the range of 0.44 to 0.46 result also in ^°Ti, 
^"^Cr, ^^Fe, and ^^'^^Ni. **Ca with Z/A»0.42 is only produced if Y^ approaches values smaller than 0.44 (see Hartmann 
et al. 1985). 
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Fig. 10. — Same as ||for the cases CS50 and W7. One sees that the low region becomes more extended with increasing 
burning front velocities, while the central values increase somewhat. ^^Co and ^^Fe are not plotted here for W7 but are 



are shown in Figure 14 
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Fig. 11. — Same as @ for the cases WSL and WLAM. WLAM has an almost constant Ye{Mr) with values close to those 
of WSL at 0.05 M© (see also Figure |). 




Mass Number A 



Fig. 12. — Ratio of abundances to solar predicted in model W7 (this is a recalculation of the 1986 model [Thielemann, 
Nomoto & Yokoi 1986] with presently available updated reaction rates and a screened NSE treatment for temperatures 
beyond 6 x lO^K, as described in Hix & Thielemann 1996). Isotopes of one element are connected by lines. The ordinate is 
normalized to ^^Fe. Intermediate mass elements exist, but are underproduced by a factor of 2-3 in comparison to Fe-group 
elements. Among the Fe-group, ^^Cr and ^^Ni are overproduced by a factor of 4, which exceeds the permitted factor of 
~ 3. 




Fig. 13. — Same as Figure |T2| i.e. a deflagration model treated exactly like W7, but with an initial composition 
corresponding to zero metallicity (50% mass fractions of ^^C and ^^O without ^^Ne admixture, i.e. Ye— 0.5). The 
production of neutron-rich species ^'*Fe, ^^Ni, and ^^Ni are smaller than the original W7 model. Noticeable is also the 
almost complete absence of intermediate odd-Z elements which are produced in explosive oxygen burning for a Ye smaller 
than 0.5. The even-Z intermediate mass elements are also effected. They are not set off by a constant factor in comparison 
to ^^Fe, as seen in Figure 12 for W7, but increase with mass. The increasing deviation between N = Z and a line in the 
nuclear chart with Ye < 0.5, with increasing nuclear mass, does not occur. 
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Fig. 14. — Composition of W7 and WTO. The major changes are given in the mass zones outside the central O.SMq, 

where Ye is cnhcritcd from the initial mctallicity in form of ^^Nc and not due to electron captures. This has effects on the 
Fe-group composition in the alpha-rich freeze-out region (on •'^^'•'^'''Ni and '^^Fe, the latter two decaying to ^^Fe and ^^Cr) 
and the incomplete Si-burning region (^^'^^Ni, ^^Co, and ^^'^^Fe, the unstable species decaying to '^^Fc, "'''''Mn, and ^^Cr). 
58,57]\jj^ ^^Co, and ^''Fe decrease below the plot limits in the incomplete burning region while ^^Fe (N=Z) increases. 
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Fig. 15. — Comparison to solar for SN la compositions which arc made up of the slow deflagrations CS15 and WS15 in 
the central layers and a W7 composition in the outer layers. This is mainly a test for the Fe-group composition, especially 
SO-pj^ 54Cj., 58Fg^ and ***Ni. 
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Fig. 16.-- Comparison to solar for a SN la composition made up from CS30 and WS30 in the central layers and W7 in 
the outer layers. The neutron-rich species decline in comparison with CS15 and WS15. 
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Fig. 17. — Comparison to solar for CS50 with a W7 composition in the outer layers. A further decrease of (very) neutron- 
rich species is visible, but ^^Fe and '^^Ni, originating from intermediate Ye regions, increase to excessive values due to the 
flatter Ye-gradient of faster deflagrations (see Figure ||). 
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Fig. 18. — Comparison to solar for WSL and WLAM with a W7 composition in the outer layers. Excessive ratios for 
neutron-rich species are observed due to extended low Ye-regions out to 0.06 Mq (see Figure ^). 
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Fig. 19. — Composition of WS15DD1 and WS15DD2 against the expansion velocity and M(r). 
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Fig. 20.— Composition of WS15DD3. In the scries DD1-DD3 we sec a decrease in the total amount of intermediate mass 
elements (Si-Ca), an increase in ^^Ni, and a change of the ratio between matter experiencing an alpha-rich freeze-out 
(indicated by the ^^Ni-plateau) and incomplete Si- burning (^^Fe-plateau). 
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Fig. 21 — Composition of CS15DD1 and CS15DD2. The outer layers affected by the detonation behave almost identical 
to the WSDD series. 
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Fig. 22. — Comparison to solar for CS15DD1 and CS15DD2. A decrease of the Si-Ca/Fe ratio is observable, a decrease 
of e.g. ^^Mn and ^^Cr (dominated by incomplete Si-burning), and an increase of ^^Ni (decaying from ^^Zn) and ^^Co 
(alpha-rich freeze-out). 
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Fig. 23.— Comparison to solar for WS15DD1 and WS15DD2. 
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Fig. 24. — Comparison to solar for WS15DD3. In this series DD1-DD3 we see similar changes as shown in Figure P4 
Differences in the neutron-rich species ^°Ti, ^''Cr, ^^Fe, and ^^Ni are due to the differences between the W and C series 
in the central slow deflagration layers (higher ignition densities). 
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Fig. 25. - Mass fractions in the models CS15DD1-DD3 as a function of expansion velocity, which reveals more easily 
the outer layers affected by the deflagration-detonation transition and the quenching of nuclear burning up to unburned 
matter in the very outer layers. Intermediate mass elements beyond Mg are essentially unburned for velocities larger 
than 20000-30000 km and would be found in their solar values scaled with the appropriate metallcity. The initial 
composition of the white dwarf for the present calculation consists only of ^^C, ^^O, and ^^Ne (according to the CNO 
metallicity of the accreted matter). 
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Table 1: Masses of ^^Ni in SN la Models Investigated 



Model 


Mass 




Energy 


''"Ni Mass 


W7 


1.38 Mq 


1.3x 


lO'^i 


ergs 


0.59 Mq 


W70 


1.38 Mq 


1.3x 


10^1 


ergs 


0.64 Mq 


WS15DD1 


1.38 Mq 


1.33X 


10^1 


ergs 


0.56 Mq 


WS15DD2 


1.38 Mq 


1.40X 


1051 


ergs 


0.69 Mq 


WS15DD3 


1.38 M0 


1.43X 


1051 


ergs 


0.77 Mq 


CS15DD1 


1.38 M0 


1.34X 


1051 


ergs 


0.56 Mq 


CS15DD2 


1.38 M0 


1.44 X 


1051 


ergs 


0.74 Mq 



Table 2: Nucleosynthesis of Neutron-rich Species (Xi/X(5^Fe)/Solar Ratio) 



Model 




P9 


Vdef/Vs 


Y 

^ c.c 






b«Fe 




"^Ni 


54Fc 


=^«Ni 


WSL 


2, 


,12 


~0.01 


0.442 


39 


74 


15 


7 


5 


1.2 


2.1 


WLAM 


2, 


,12 


~0.01 


0.447 


22 


47 


9 


1.6 


5 


1.8 


3.8 


WS15W7 


2, 


.12 


0.015 


0.440 


4 


7 


1.7 


0.7 


2.1 


1.3 


2.7 


WS15DD1 


2, 


.12 


0.015 


0.440 


4 


7 


1.8 


0.7 


1.1 


1.9 


1.3 


WS15DD2 


2, 


,12 


0.015 


0.440 


3.1 


6 


1.5 


0.6 


2.4 


1.3 


1.5 


WS15DD3 


2, 


,12 


0.015 


0.440 


2.8 


5 


1.3 


0.5 


3.2 


1.0 


1.5 


WS30W7 


2, 


,12 


0.03 


0.445 


1.3 


4 


1.5 


0.1 


2.5 


1.3 


2.6 


CS15W7 


1, 


.37 


0.015 


0.449 


0.5 


2 


0.9 


0.01 


1.9 


1.3 


2.7 


CS15DD1 


1, 


.37 


0.015 


0.449 


0.2 


1.1 


0.5 


0.01 


0.5 


2.0 


1.4 


CS15DD2 


1, 


,37 


0.015 


0.449 


0.6 


2.7 


1.3 


0.01 


2.8 


1.2 


1.4 


CS30W7 


1, 


,37 


0.03 


0.455 


0.07 


0.5 


0.3 


0.002 


2.0 


1.4 


2.7 


CS50W7 


1, 


,37 


0.05 


0.459 


0.003 


0.04 


0.04 


0.0004 


1.8 


1.8 


3.2 


W7 


2, 


.12 


up to 0.3 


0.446 


1.1 


3.8 


1.6 


0.05 


2.3 


2.4 


4.3 


W70 


2, 


.12 


up to 0.3 


0.446 


0.97 


3.4 


1.4 


0.04 


1.3 


1.7 


3.0 
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Table 3: Nucleosynthesis products of SN II and la models 







Synthesized mass 


(Mo) 






lype 11 
Species lO-SOM© 


W7 


WTO WDDl 


Type la 

WDD2 WDD3 


CDDl 


CDD2 





7, 


,93E-02 


4 


.83E- 


-02 


5, 


.08E- 


-02 


5, 


,42E- 


-03 


8, 


,99E- 


-03 


1 


.66E- 


-02 


9, 


.93E- 


-03 


5, 


.08E- 


-03 


13c 


3 


.80E-09 


1 


.40E- 


-06 


1, 


.56E- 


-09 


5, 


.06E- 


-07 


3 


.30E- 


-07 


3 


.17E- 


-07 


8, 


.46E- 


-07 


4, 


.16E- 


-07 




1 


.56E-03 


1 


.16E- 


■06 


3. 


.31E- 


■08 


2, 


.84E- 


■04 


2, 


.69E- 


•04 


1 


.82E- 


•04 


9, 


.06E- 


•05 


9, 


.03E- 


■05 


15N 


1 


.66E-08 


1 


.32E- 


■09 


4, 


.13E- 


■07 


9, 


.99E- 


■07 


5, 


.32E- 


•07 


1 


.21E- 


•07 


2, 


.53E- 


•07 


2, 


.47E- 


■07 


160 




1.80 


1 


.43E- 


-01 


1, 


.33E- 


-01 


8, 


.82E- 


-02 


6 


.58E- 


-02 


5, 


.58E- 


-02 


9, 


.34E- 


-02 


5, 


.83E- 


■02 




9 


.88E-08 


3 


.54E- 


■08 


3, 


.33E- 


•10 


3, 


.77E- 


•06 


4, 


.58E- 


•06 


3 


.60E- 


•06 


9, 


.55E- 


•07 


1, 


.OlE- 


•06 


18Q 


4, 


.61E-03 


8, 


.25E- 


■10 


2, 


.69E- 


■10 


6, 


.88E- 


■07 


6, 


.35E- 


•07 


2 


.39E- 


•07 


2, 


.08E- 


•07 


1, 


.92E- 


•07 


19p 


1 


.16E-09 


5 


.67E- 


-10 


1, 


.37E- 


■10 


1, 


.70E- 


■09 


4, 


.50E- 


■10 


2 


.30E- 


-10 


5, 


.83E- 


-10 


4, 


.24E- 


-10 


20Ne 


2 


.12E-01 


2, 


.02E- 


-03 


2, 


,29E- 


-03 


1, 


.29E- 


-03 


6, 


.22E- 


-04 


4, 


.55E- 


-04 


1, 


,16E- 


-03 


6, 


.05E- 


-04 


2iNe 


1 


.08E-03 


8, 


.46E- 


■06 


2 


.81E- 


■08 


1, 


.16E- 


■05 


1, 


.39E- 


■06 


1 


.72E- 


•06 


3. 


.63E- 


■06 


1, 


.99E- 


■06 


22Ne 


1 


,83E-02 


2, 


,49E- 


-03 


2 


.15E- 


■08 


1, 


.51E- 


■04 


4, 


.21E- 


■04 


8 


.25E- 


-04 


4, 


,41E- 


-04 


2 


.llE- 


■04 


23Na 


6 


,51E-03 


6, 


.32E- 


-05 


1, 


.41E- 


■05 


8, 


.77E- 


■05 


2, 


.61E- 


■05 


3 


.OlE- 


-05 


5, 


lOE- 


-05 


3, 


.50E- 


■05 


24Mg 


8, 


.83E-02 


8, 


.50E- 


-03 


1, 


.58E- 


■02 


7, 


.55E- 


■03 


4, 


.47E- 


-03 


2, 


.62E- 


-03 


7, 


,72E- 


-03 


4, 


,20E- 


-03 


25Mg 


1 


.44E-02 


4, 


.05E- 


■05 


1, 


.64E- 


■07 


8, 


.23E- 


•05 


2, 


.66E- 


■05 


2 


.68E- 


•05 


4, 


.85E- 


■05 


3, 


.25E- 


■05 


26Mg 


2 


.OlE-02 


3 


.18E- 


■05 


1, 


.87E- 


■07 


6, 


.25E- 


•05 


2, 


.59E- 


■05 


1 


.41E- 


•05 


4, 


,96E- 


■05 


2, 


.97E- 


■05 


27 Al 


1 


.48E-02 


9 


.86E- 


-04 


1, 


.13E- 


■04 


4, 


.38E- 


■04 


2, 


,47E- 


-04 


1 


.41E- 


-04 


4, 


,45E- 


-04 


2, 


.35E- 


■04 


28Si 


1 


.05E-01 


1 


.54E- 


■01 


1, 


.42E- 


-01 


2, 


.72E- 


-01 


2, 


.06E- 


•01 


1 


.58E- 


•01 


2, 


,77E- 


■01 


1, 


.98E- 


-01 


29Si 


8, 


.99E-03 


9 


.08E- 


■04 


5, 


.79E- 


■05 


5, 


.47E- 


•04 


3 


.40E- 


•04 


2 


.13E- 


•04 


5, 


.52E- 


■04 


3 


.22E- 


■04 


30Si 


8 


.05E-03 


1 


.69E- 


-03 


7, 


.12E- 


■05 


1, 


.03E- 


■03 


6, 


.41E- 


■04 


3 


.88E- 


-04 


1, 


,05E- 


-03 


6 


,14E- 


-04 


31p 


1 


.21E-03 


3 


.57E- 


-04 


9, 


,12E- 


-05 


2, 


.38E- 


■04 


1, 


.60E- 


-04 


1 


.07E- 


-04 


2, 


.40E- 


-04 


1, 


,52E- 


-04 


32g 


3 


.84E-02 


8, 


.46E- 


■02 


9, 


.14E- 


•02 


1, 


.60E- 


-01 


1, 


,22E- 


-01 


9 


.37E- 


-02 


1, 


.63E- 


-01 


1, 


,17E- 


-01 


33s 


1 


.78E-04 


4 


,24E- 


-04 


6 


.07E- 


■05 


2, 


.74E- 


■04 


1, 


,92E- 


-04 


1 


.34E- 


-04 


2, 


.71E- 


-04 


1, 


.79E- 


■04 


34g 


2 


.62E-03 


1 


.98E- 


-03 


1, 


.74E- 


■05 


2, 


.76E- 


■03 


2, 


,04E- 


-03 


1 


.46E- 


-03 


2, 


.77E- 


-03 


1, 


.90E- 


■03 


36 g 


1, 


.78E-06 


4, 


.18E- 


-07 


3 


.41E- 


■11 


2, 


.23E- 


-07 


1, 


,31E- 


-07 


7, 


.44E- 


-08 


2, 


.22E- 


-07 


1, 


,25E- 


■07 


35C1 


1 


.OlE-04 


1 


.37E- 


■04 


1, 


.06E- 


•05 


9, 


.28E- 


•05 


7, 


.07E- 


•05 


5, 


.33E- 


•05 


9, 


.03E- 


■05 


6, 


.56E- 


■05 




1 


.88E-05 


3 


.67E- 


■05 


5, 


.56E- 


■06 


2, 


.94E- 


•05 


2, 


.26E- 


•05 


1 


.71E- 


•05 


2, 


,86E- 


■05 


2, 


llE- 


■05 


36 Ar 


6 


.62E-03 


1 


.47E- 


-02 


1, 


,91E- 


-02 


3 


.llE- 


■02 


2, 


,41E- 


-02 


1 


.87E- 


-02 


3, 


,18E- 


-02 


2, 


.34E- 


-02 


38 Ar 


1 


.37E-03 


9 


.50E- 


-04 


6, 


.60E- 


■07 


1, 


.23E- 


■03 


9, 


.90E- 


■04 


7 


.44E- 


-04 


1, 


,20E- 


-03 


9, 


.26E- 


■04 


4"Ar 


2 


.27E-08 


1 


.87E- 


-08 


3, 


.42E- 


■12 


1 


.81E- 


■09 


•J 


.19E- 


■09 


3 


..jGE- 


-09 


/ , 


.•■)GE- 


-09 


4, 


.82E- 


■09 
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Continued form TABLE 3. 



Synthesized mass [Mq) 



Species 


'J 

10 


-ype 11 
i-5OM0 




W7 






W70 


WDDl 


'I'ype la 
WDD2 


WDD3 


CDDl 


CDD2 




6, 


.23E- 


-05 


7, 


.23E- 


■05 


1, 


.67E- 


-06 


6 


.76E- 


-05 


5, 


.67E-05 


4, 


.52E- 


-05 


6 


.39E- 


-05 


5, 


,34E- 


-05 




5, 


.07E- 


■06 


6 


.llE- 


•06 


4, 


.83E- 


•07 


5, 


.43E- 


•06 


4, 


.52E-06 


3, 


.62E- 


•06 


5, 


.20E- 


•06 


4, 


.25E- 


•06 


40Ca 


5, 


,77E- 


■03 


1 


.19E- 


■02 


1, 


.81E- 


■02 


3 


lOE- 


■02 


2, 


.43E-02 


1, 


.88E- 


■02 


3 


.18E- 


■02 


2, 


.38E- 


■02 


42Ca 


4, 


,23E- 


-05 


2 


.82E- 


■05 


1, 


.06E- 


■08 


3 


.09E- 


■05 


2, 


.55E-05 


1, 


.93E- 


-05 


2, 


.97E- 


-05 


2, 


.36E- 


■05 


43Ca 


1, 


,08E- 


•06 


9 


.64E- 


■08 


6, 


.17E- 


-08 


6, 


.60E- 


-08 


2, 


.22E-07 


4. 


.18E- 


-07 


5, 


.15E- 


-08 


2, 


,96E- 


-07 


44Ca 


5. 


.53E- 


■05 


8, 


.02E- 


■06 


1, 


.38E- 


■05 


1, 


.44E- 


■05 


2. 


.95E-05 


4, 


,66E- 


■05 


1 


.37E- 


■05 


3 


.62E- 


■05 


^^Ca 


1, 


,43E- 


-10 


4 


.16E- 


■09 


1, 


,01E- 


-09 


5, 


,01E- 


-09 


4, 


.73E-09 


4, 


.47E- 


-09 


8, 


.79E- 


-10 


1, 


,18E- 


-09 


^^Ca 


5, 


,33E- 


-14 


2 


.60E- 


■09 


2, 


,47E- 


-09 


1, 


.63E- 


-06 


1, 


,64E-06 


1, 


.55E- 


-06 


3 


,54E- 


-11 


4, 


.93E- 


-10 


45Sc 


2, 


,29E- 


-07 


2 


.21E- 


■07 


3 


.85E- 


■08 


2, 


.49E- 


■07 


2, 


.09E-07 


1, 


.76E- 


-07 


2, 


,47E- 


-07 


2, 


.02E- 


■07 


46Ti 


7. 


.48E- 


■06 


1 


.33E- 


■05 


3, 


.49E- 


■07 


1, 


.34E- 


■05 


1, 


.12E-05 


8, 


.58E- 


■06 


1 


.27E- 


■05 


1, 


.05E- 


■05 


47Ti 


2, 


llE- 


■06 


5 


.lOE- 


■07 


4, 


.08E- 


■07 


5, 


.65E- 


■07 


1, 


.56E-06 


2, 


.57E- 


■06 


4, 


.93E- 


■07 


1, 


.95E- 


■06 


48Ti 


1, 


,16E- 


-04 


2 


.05E- 


■04 


3 


,13E- 


■04 


7 


.lOE- 


■04 


6, 


.llE-04 


5, 


.23E- 


-04 


7 


,32E- 


-04 


6, 


.20E- 


■04 


49Ti 


5, 


,98E- 


•06 


1 


.71E- 


•05 


2, 


.94E- 


•06 


5, 


.27E- 


•05 


4, 


.39E-05 


3, 


.59E- 


•05 


5, 


.22E- 


•05 


4, 


.17E- 


•05 




3, 


,81E- 


■10 


1 


.07E- 


■04 


1, 


.04E- 


■04 


3, 


.52E- 


■04 


3 


.51E-04 


3, 


.51E- 


■04 


2 


.08E- 


■05 


7, 


.28E- 


■05 


50v 


7, 


,25E- 


-10 


1 


.55E- 


■08 


1, 


.22E- 


■08 


9, 


.74E- 


■09 


9, 


.33E-09 


9, 


.07E- 


-09 


4, 


.94E- 


-09 


1, 


,20E- 


-08 


Sly 


1, 


,00E- 


-05 


7 


.49E- 


■05 


4, 


,27E- 


-05 


1, 


.33E- 


-04 


1, 


,16E-04 


1, 


.02E- 


-04 


1 


.llE- 


•04 


1, 


.09E- 


-04 


50Cr 


4, 


,64E- 


■05 


2 


.73E- 


■04 


6, 


.65E- 


■05 


4, 


.44E- 


■04 


3 


.53E-04 


2, 


.84E- 


■04 


4, 


.49E- 


■04 


3 


.36E- 


■04 


52Cr 


1, 


.15E- 


-03 


6 


.36E- 


03 


7, 


.73E- 


-03 


1, 


.68E- 


■02 


1, 


.37E-02 


1, 


,13E- 


-02 


1, 


.65E- 


-02 


1, 


.40E- 


-02 


53 Cr 


1, 


.19E- 


-04 


9 


.22E- 


■04 


5, 


.66E- 


■04 


1, 


.66E- 


■03 


1, 


.38E-03 


1, 


.17E- 


-03 


1, 


.59E- 


-03 


1, 


.33E- 


-03 


54Cr 


2, 


.33E- 


-08 


9 


.24E- 


■04 


9, 


.04E- 


■04 


1, 


.60E- 


■03 


1, 


.60E-03 


1, 


.60E- 


-03 


2, 


,31E- 


-04 


8, 


,02E- 


-04 


55Mn 


3, 


.86E- 


•04 


8, 


.87E- 


•03 


6, 


.66E- 


■03 


8, 


.48E- 


■03 


7, 


.05E-03 


6, 


.16E- 


•03 


8, 


.lOE- 


■03 


6, 


.77E- 


■03 


54Fe 


3, 


,62E- 


•03 


9 


.55E- 


■02 


7, 


.30E- 


■02 


7, 


.08E- 


■02 


5, 


,91E-02 


5, 


.15E- 


■02 


7, 


.20E- 


■02 


5, 


.64E- 


■02 




8, 


,44E- 


-02 


6 


.26E- 


■01 


6, 


.80E- 


■01 


5, 


.87E- 


■01 


7 


.13E-01 


7, 


,95E- 


-01 


5, 


.65E- 


-01 


7 


.57E- 


-01 


57Fe 


2, 


,72E- 


-03 


2 


.45E- 


•02 


1, 


.92E- 


-02 


1, 


.08E- 


-02 


1, 


.67E-02 


2, 


.06E- 


-02 


1 


.OlE- 


•02 


1, 


.80E- 


•02 




7, 


.22E- 


■09 


3 


.03E- 


■03 


2, 


.96E- 


■03 


3 


.23E- 


■03 


3 


.23E-03 


3, 


.24E- 


■03 


8, 


.63E- 


■04 


3 


.06E- 


■03 


59C0 


7, 


,27E- 


-05 


1 


.04E- 


03 


9, 


.68E- 


■04 


3 


.95E- 


-04 


6 


,25E-04 


7, 


.75E- 


-04 


2, 


,91E- 


-04 


6 


.35E- 


-04 


58Ni 


3, 


.63E- 


-03 


1, 


.lOE- 


■01 


8, 


,34E- 


-02 


3 


,14E- 


-02 


4, 


,29E-02 


4, 


.97E- 


-02 


3 


.15E- 


-02 


4, 


,47E- 


-02 


6"Ni 


1, 


.75E- 


■03 


1 


.24E- 


02 


1, 


,47E- 


-02 


5, 


.08E- 


-03 


1, 


,15E-02 


1, 


.67E- 


•02 


2 


.81E- 


-03 


1, 


,21E- 


-02 


61 Ni 


8, 


.35E- 


■05 


2 


.35E- 


■04 


2, 


.15E- 


■04 


7, 


.OOE- 


■05 


3, 


.58E-04 


5, 


.92E- 


■04 


4, 


.OOE- 


■05 


4, 


.27E- 


■04 



Continued form TABLE 3. 









Synthesized mass 


(Mq) 








Species 


Type II 
1O-5OM0 


W7 


W70 


WDDl 


Type la 
WDD2 


WDD3 


CDDl 


CDD2 


a^Ni 


5.09E-04 


3.07E-03 


1.85E-03 


1.37E-03 


3.69E-03 


5.46E-03 


6.51E-04 


4.60E-03 


64Ni 


3.20E-14 


1.70E-05 


1.65E-05 


2.32E-04 


2.31E-04 


2.32E-04 


2.47E-06 


9.29E-06 


63CU 


8.37E-07 


2.32E-06 


3.00E-06 


2.77E-06 


4.88E-06 


5.92E-06 


6.14E-07 


3.97E-06 


65CU 


4.07E-07 


6.84E-07 


8.33E-07 


7.08E-07 


2.04E-06 


3.38E-06 


1.14E-07 


2.05E-06 


64Zn 


1.03E-05 


1.06E-05 


7.01E-05 


3.71E-06 


3.10E-05 


5.76E-05 


1.87E-06 


3.96E-05 


66Zn 


8.61E-06 


1.76E-05 


6.26E-06 


2.16E-05 


6.42E-05 


1.04E-04 


2.84E-06 


6.11E-05 


67Zn 


1.52E-08 


1.58E-08 


7.28E-09 


6.35E-07 


6.55E-07 


6.27E-07 


1.69E-09 


4.01E-08 


68Zn 


3.92E-09 


1.74E-08 


1.13E-08 


7.44E-08 


8.81E-08 


9.42E-08 


3.08E-09 


3.03E-08 
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Table 4: Major Radioactive Species in SN la Models 



Syiilliosiz(Hl masri (M,., ) 


Species 


W7 


WTO 


Type la 
WDDl WDD2 WDD3 


CDDl 


CDD2 



IN a 


-1 

i. 


.YorLr 


no 

-Oo 


i 


.vobj- 


-Uo 


o 
/ 


.OOrj- 


-U» 


-1 
i 


.66tj- 


-Oo 


o 
3, 


.U4rj- 


AO 

-Oo 


i 


^? OT7 


-Uo 


y. 


64E-09 


26 Al 


4. 


.93E- 


■OT 


2, 


.92E- 


•08 


4, 


.61E- 


•OT 


1 


.88E- 


•OT 


1, 


.64E- 


■OT 


3. 


.81E- 


■OT 


2, 


.32E- 


■OT 


36Q 


2, 


.58E- 


•06 


3 


.9TE- 


•10 


1, 


.25E- 


•06 


7 


.85E- 


•OT 


4, 


.89E- 


■OT 


1, 


.23E- 


■06 


T, 


.35E- 


■OT 


39 Ar 


1, 


.20E- 


■08 


2, 


•OOE- 


-13 


6, 


.04E- 


-09 


4 


.44E- 


-09 


3, 


.29E- 


-09 


5, 


.99E- 


-09 


4, 


.12E- 


-09 


40k 


8, 


.44E- 


-08 


5, 


.46E- 


•12 


3, 


.29E- 


-08 


2 


.35E- 


■08 


1, 


.T3E- 


■08 


3, 


.13E- 


-08 


2, 


.15E- 


■08 


4iCa 


6, 


,09E- 


•06 


4, 


.83E- 


•OT 


5, 


.42E- 


•06 


4, 


.52E- 


•06 


3, 


.62E- 


■06 


5, 


.20E- 


■06 


4, 


.25E- 


■06 




T, 


.94E- 


■06 


1, 


.38E- 


-05 


1, 


.44E- 


-05 


2 


.95E- 


-05 


4, 


.65E- 


-05 


1, 


.36E- 


-05 


3, 


.62E- 


-05 


48v 


4, 


.95E- 


-08 


1, 


.61E- 


-08 


6, 


.52E- 


-08 


5, 


.38E- 


■08 


4, 


.llE- 


■08 


6, 


.33E- 


-08 


5, 


.03E- 


-08 


49v 


1, 


.52E- 


■07 


3, 


.23E- 


•08 


1, 


.18E- 


•OT 


1 


.03E- 


■OT 


8, 


.63E- 


■08 


9, 


.99E- 


■08 


8, 


.88E- 


■08 


53Mn 


2, 


.TTE- 


■04 


2, 


.48E- 


-04 


1, 


.58E- 


-04 


1 


.54E- 


-04 


1, 


.50E- 


-04 


8, 


.4TE- 


-05 


9, 


.69E- 


-05 


60 pg 


T, 


.52E- 


-OT 


T, 


.19E- 


-OT 


7 


.33E- 


-05 


7 


.33E- 


-05 


T, 


.36E- 


-05 


6, 


.23E- 


-08 


2, 


.T4E- 


-OT 


56Co 


1, 


.44E- 


-04 


1, 


,26E- 


-04 


6 


,04E- 


-05 


5 


,51E- 


-05 


5, 


.20E- 


-05 


6, 


.18E- 


-05 


5, 


,24E- 


-05 


s^Co 


1, 


.48E- 


•03 


1, 


.42E- 


•03 


6, 


.43E- 


•04 


6 


.42E- 


■04 


6, 


.36E- 


■04 


4, 


.50E- 


•04 


4, 


.4TE- 


■04 


60Co 


4, 


,22E- 


•OT 


4, 


.13E- 


•OT 


3 


.T6E- 


•OT 


3 


.81E- 


■OT 


3, 


.T8E- 


■OT 


1, 


.09E- 


■OT 


4, 


.23E^ 


■OT 


56Ni 


5, 


.86E- 


-01 


6, 


.41E- 


-01 


5, 


.64E- 


-01 


6 


.90E- 


-01 


T, 


.T3E- 


-01 


5, 


.55E- 


-01 


T, 


.3TE- 


-01 


"Ni 


2, 


.2TE- 


-02 


1, 


.T5E- 


-02 


9, 


.95E- 


-03 


1 


.59E- 


-02 


1, 


.98E- 


■02 


9, 


.5TE- 


•03 


1, 


.T4E- 


■02 


59Ni 


6, 


.TIE- 


•04 


6, 


.39E- 


•04 


2, 


.54E- 


•04 


2 


.53E- 


■04 


2, 


.51E- 


■04 


2, 


.03E- 


■04 


1, 


.90E- 


■04 


63Ni 


8. 


.OOE- 


•OT 


T. 


.81E- 


•OT 


1. 


.69E- 


•06 


1 


.T4E- 


■06 


1, 


.T3E- 


■06 


1, 


.98E- 


■OT 


8, 


,31E- 


■OT 



Table 5: N^/Nu Ratios obtained for different Models of SN la 



SN II Model 


WT 


WTO 


DDI 


DD2 


DD3 


20 Mq 
10-50 Mq 


0.214 
0.295 


0.1T6 
0.243 


1.08 
1.48 


0.18T 
0.25T 


0.146 
0.201 



